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Abstract 

Semiclassical wave functions in billiards based on the Maslov - Fedoriuk approach 
are constructed. They are defined on classical constructions called skeletons which are 
the billiards generalization of Arnold's tori. Skeletons in the rational polygon billiards 
considered in the phase space can be closed with a definite genus or can be open being 
a cylinder-like or Mobius-like bands. The skeleton formulation is applied to calculate 
semiclassical wave functions and the corresponding energy spectra for the integrable and 
pseudointegrable billiards as well as in the limiting forms in some cases of chaotic ones. 
The superscars of Bogomolny and Schmit are shown to be simply singular semiclassical 
solutions of the eigenvalue problem in the billiards well built on the singular skeletons in 
the billiards with fiat boundaries in both the integrable and the pseudointegrable billiards 
as well as in the chaotic cases of such billiards. 

PACS number(s): 03.65.-w, 03.65.Sq, 02.30.Jr, 02.30.Lt, 02.30.Mv 

Key Words: Schrddinger equation, semiclassical expansion, Lagrange manifolds, classical 
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1 Introduction 



Billiards while a non-analytic motion area are however well known as examples of the 
non-integrable two dimensional systems except the known cases of the integrable elliptical, 
rectangular and some triangle billiards. They are widely considered as a simple field of exper- 
imental [1, 12, 3] as well as theoretical [4, 5, 6, 7] (and papers cited there) and computational 
investigations [8, 9, 10] allowing to apply many different methods (see Sarnak's lecture [11] 
and [12] of the same author for an extensive review of the respective theoretical methods 
covering also billiards manifolds as well as the students book of Tabachnikov [13]). 

In this paper we are going to develop the semiclassical wave function (SWF) formalism 
which can be applied to non-integrable cases of the two dimensional motions in billiards and 
which can be easily extended to higher dimensions. 

Essentially our approach is very close to the one of Maslov and Fedoriuk [14]. The main 
difference between them is in a treatment of crossing the singular points of the SWF's set 
on caustics. Namely, instead of making the canonical phase space variable transformations 
accompanied by the Fourier transformations of the SWF's to move through the caustic points 
we apply the analytical continuation on the complex time plane to both the SWF's and the 
classical trajectories. This greatly simplifies the corresponding procedure in comparison with 
the Maslov and Fedoriuk treatment and allows us for not leaving the configuration space [22] . 
However in this paper there is no opportunity to use the simplification mentioned. 

It is further a classical construction which we called skeletons on which the SWF's are 
defined. Each skeleton is compound in a closed way of bundles of rays (classical trajectories). 
SWF's of basic forms are defined just first on bundles while the global SWF's (GSWF) 
satisfying vanishing boundary conditions of Dirichlet or Neumann are sums of these basic 
SWF's (BSWF) and are uniquely and continuously defined on the skeleton. Skeletons play 
in this way a role of Arnold's tori [15] except that a number of ray bundles in skeletons can 
be infinite if billiards with chaotic motions are considered. 

When forms of billiards boundaries are considered and their relationships with types of 
the classical motion in them one realizes that billiards can be dividing into two general classes: 
the one with the flat boundaries, i.e. the polygon billiards and the second class in which some 
pieces of their boundaries have finite curvature. In the polygon billiards there are no caustic 
phenomenon and this is the main and essential property which differs both the classes. 

However among the polygon billiards one can still distinguish a class of the rational bil- 
liards, i.e. which all angles are rational part of vr. As it was shown by Richens and Berry 
[8] (see also Tabachnikov [13]) the rational billiards can be classified according to classical 
motions in the corresponding phase space which is performed on a set of two-dimensional 
disjoint compact Lagrangian surfaces which collecting together can be made equivalent topo- 
logically to a compact close surface of a given genus defined by the billiards. All possible 
values of genus can be realized by the rational billiards. 

As an archetype of the rational pseudointegrable polygon billiards can be considered the 
broken rectangular billiards, i.e. billiards which can be glued of a finite number of rectangles. 
The SWF's in the rectangular billiards and in the broken ones are just studied intensively 
in this paper. Nevertheless other simple polygons such as the equilateral triangles and the 
pentagon billiards arc also considered. 

It is well known that the semiclassical limits of wave functions are typically singular, i.e. 
the semiclassical wave functions are frequently deprived of such properties as finiteness and 
smoothness which have to be satisfied by the exact wave functions. Nevertheless SWF's alone 
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as well as accompanied them energy spectra give typically very good approximations to the 
exact ones in the high energy limit. 

However it happens frequently that semiclassical calculations provide us with exact wave 
functions and energy spectra. This can takes place for example when the calculated SWF's 
satisfy the same conditions as the exact ones and the classically allowed configuration space 
coincides with the quantum one. Such SWF's are called regular in this paper. Examples 
of such regular semiclassical solutions are provided by the generic SWF in the rectangle, the 
equilateral triangle and the pentagon billiards. 

Most of the semiclassical wave functions which solve the eigenvalue problem in the semi- 
classical limit are however singular, i.e. they do not satisfy some of the conditions which 
the exact solutions have to do. The SWF's which can be constructed for the rectangular 
billiards and the broken rectangle ones arc not exceptional in the respects mentioned and can 
be also of two types - the regular and the singular ones depending on skeletons used to their 
constructions. It is just the singular SWF's for which superscars phenomenon of Bogomolny 
et al [17, 18] can be observed. While this singular behaviour of SWF's in the cases considered 
is reduced only to discontinuities of the first derivatives of the respective SWF's this "defect" 
of them is enough for being a source of the superscar phenomenon. 

It is also shown in the paper that the superscars are typical not only for the integrable 
and pseudointegrable systems but also for the chaotic ones. A well known example is the 
Bunimovich stadium with its bouncing ball modes [3, 4, 5]. But it is easy to give many other 
examples of the chaotic billiards with any form of the superscars which can be found in the 
broken rectangular billiards. 

The paper is organized as follows. 

In the next section the Maslov - Fedoriuk method of the semiclassical wave function 
construction is reminded and discussed. 

In sec. 3 a construction of skeletons is given. 

In sec.4 global SWF's in billiards are constructed. 

In sec. 5 the rectangular billiards, the broken rectangle ones, the triangle and the pentagon 
billiards are considered. In the rectangular billiards case all possible SWF's which can be 
defined in it, i.e. the regular and the singular ones, are discussed. SWF's in the remaining 
billiards are considered selectively because of increasing complexity of the corresponding 
skeletons. Nevertheless it is shown that in these cases the superscar SWF's are common and 
more spectacular than for the rectangular billiards. 

Next in sec. 6 it is shown that the superscar SWF's can be implemented into chaotic 
deformations of the broken rectangle, the triangle and the pentagon billiards. 

In sec.7 the results of the paper are summarized. 

There are two appendixes attached to the paper which justify the main assumptions used 
in the construction of the global SWF's on skeletons. 

2 Semiclassical wave function expansion for n-D stationary 
Schrodinger equation 

Consider the n-dimensional stationary Schrodinger equation: 

A^(r) + A2|j(E-y(r))*(r) =0 (1) 
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with a potential V{r), r € i?„ confining a point particle with a mass m and containing a 
formal dimensionless parameter A > 0. For a convenience we shall put further h = 1 and 
m = 1. The Schrodinger equation is recovered by putting A = 1 in (1). 

We would like to construct a solution to Eq.(l) using the idea of Maslov et al [14] and 
considering the wave function \l'(r) as defined on families of classical trajectories a dynamic 
of which is given by the classical Hamiltonian H = ^p"^ + V(r) and which carry an energy 
Eq all. 

Such families are constructed locally in the following way. 

In Rn we choose a n—l-D hypcrsurface parametrized by local coordinates (si, s„_i) 
On the initial momenta p(ro), tq G are defined so that the pair (ro,p(ro)), tq G 

serve as the initial data for the trajectories r{t) = f (ro, p(ro); i) developed by the 
Hamiltonian H. Additionally the momentum field p(ro) defined on has to satisfy: 

£p(ro)dro = (2) 

for each loop C, C C Sn-i- 

We can now define the transformation: r — t- (t, Tq) — )■ (t, si, s„_i), r= (xi, .t„), Tq = 
i), ■■■,xo,nisi, Sn-i)), which is one-to-one up to a caustic surface C„_i on 
which the Jacobean (f (ro, p(ro); i) = f (t, si, s„_i)): 



J(t,Si, ...,Sn-l) 



9^ 5^ 
9i ' dsj 



(3) 



vanishes. 

A n-dimensional domain A„ of 2ra-dimensional phase space i?2n made in this way by the 
hypcrsurface Sn-i and trajectories emerging from it is known as the Lagrange manifold [15]. 

Therefore in the variables t, si, Sn-i the new wave function si, satisfies 
the following relation with the previous one: 

|x(t,Sl,...,S„-l)P = |*(f(t,Sl,...,S„-l)|^|J(t,Sl,...,Sn-l)| (4) 

The particle momentum p on the trajectories r(t) = f (t, si, satisfies of course 
the equation: 

= p(f(t,si,...,s„_i)) (5) 

defining also the Jacobean evolution. Namely: 

d dfi{t,Si,...,Sn-l) ^ A dpi dfj{t,Si,...,Sn-l) 

dt dt ~^^dxi dt 

dpi dfj{t, 



dt dsi dxj dsi 



so that 



dJ{t,Sl, ...,Sn-l] 

dt 



/ = l,...,n-l (6) 



J{t,Si, ...,SH_i)Vp(f(t,Si, ...,Sn-l)) (7) 
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The above equation is just the Liouville theorem with the solution: 

J{t,s,,...,sn-i) = J(si,...,s„_i)e/o^P(f(*''^i.-.-"-i))'^*' (8) 

where J{si, s„_i) is the value of the Jacobean on the hypersurface Sn-i- 

It is well known from the classical Hamiltonian mechanics [15] that the action integral: 

5(r,ro)= l\{v')dv' (9) 

J TO 

taken on the Lagrange manifold A„ is a point function of r and vq. Therefore taking fq 
as a definite fixed point of the hypersurface Sn-i and denoting by 5(r) the action function 
corresponding to this case we can complete a definition of the wave function si, s„_i) 
by the following equation: 

^%i{t, si, Sn-i)) = J-^ {t, Si, s„_i))e'^^^^(f(*'^i--^"-i)x''(i, Si, Sn-l) (10) 

where cr = it is a signature of *'^(r). 

The form (10) of the semiclassical wave functions (SWF) will be called basic (BSWF). 
Therefore the quantities involved in the above definitions satisfy the following equations: 



p(r) = V5(r) 
^p2(r) + V{r) -Eo = 

A(J-^X''(r)) + (j2iAJ-^(r)Vx''(r) • p(r) + 2X^{E - Eo)J-^2{v)x'' {v) = 



r = f(t,si,...,s.„_i) (11) 

By the variables i, si, the third of the last equations can be rewritten in the 
following form: 

^^'^ dt + 

J5A(j-5x^(t,si,..,s„_i,A)) +A2(£;-Eo)x''(i,si,...,s„_i,A) = (12) 

where a dependence of x'^(^) si, A) on A was shown explicitly. 

The Eq.(12) describes the time evolution of x'^(^) ^i) ^n-i, A) along trajectories starting 
on the hypersurface Sn-i if its "initial" values on this surface, i.e. x°^(0, si, A) = 
X'^(si,---,s„_i,A) are given. 

We are going to consider the equation (12) in the semiclassical limit A — > +oo looking for 
its solutions in the form of the following asymptotic series: 



E — Eq — ^ EkX ^ ^ 
k>l 

X^it, Si, Sn-l,X) = Xkit, Si, Sn-l)\~^ 
k>0 

X^isi, sn-i, A) = ^ Xfc(si, s„_i)A-*^ (13) 

fe>0 

Putting A = 1 in (1), (10) and (13) we get approximate semiclassical solutions to the 
energy eigenvalue problem of the Schrodinger equation. 
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It is to be noticed that for the selfconsistency reasons the semiclassical series for the 
energy parameter in (13) starts from the second power of A~^, i.e. this ensures the proper 
hierarchy of steps in the algorithm of semiclassical calculations by which the higher order 
terms of the series in (13) are determined by the lower order ones. 

It should be noticed also that despite the fact that Eq enters the classical equation of 
motion (11) it is still quantum, i.e. its value depends on h which is considered to have the 
definite numerical value, i.e. H is not a parameter. In particular the series (13) represent 
the inverse power hierarchy in the formal parameter A, i.e. not in powers of H, between 
subsequent terms. 

Moreover Eq if quantized can depend on A. However, whatever this dependence is the 
semiclassical series of the difference E — Eq must be given by (13). 

Needless to say the introducing A makes a treatment of the Schrodingcr equation equiv- 
alent of course to considering it in the limit ?i — t- 0, i.e. semiclassically, clearly however 
separating the role of 7i as a parameter from its role defining the microscale of quantum 
phenomena. 

Substituting (13) into (12) we get: 

dXo{t,si,...,Sn~i) ^ „ 

dt 
at 

^ (^J^A {j-^xUt^ Si, Sn-l)) +2J2 Ek-l+lXl{t, Si, Sn-l)^ 

k = 0,1,2,..., (14) 



with the obvious solutions: 



Xo{t, Si, Sn-l) = Xoisi, Sn-l) 
Xk+l{t,Si,...,Sn-l) = Xfc+l(si,-,Sn-l) + 



y ^ i^J-2A (j-hUi', si> Sn-l)) + 2 Ek-i+ixi{t', Si, s„_i) j dt' 

k = 0,1,2,..., (15) 



For future applications it is worth to note that if (12) is obviously invariant on a reparametri- 
zation of the hypersurface Sn-i it is also invariant on the following change of variables: 

t t(si, ...,Sn-l) ± t 
Sk hkisi,...,Sn-l) 

k = l,...,n-l (16) 
if it is accompanied simultaneously by the transformations: 

X''{t,Si,...,Sn-l,X) 

X'^'^it, Si, Sn-l, A) = (±Ji)-2 (/i-l(si, Sn-l), Ki-li''^l^ «n-l)) X 

X'^(t(,Si, ,S„_i) ± t, hi^{si, Sn-l), /in-l(si> «n-l), A) (17) 

where Ji{si, Sn-i) is the Jacobean of the transformation Sk — >■ hk{si, Sn-i), k = 
l,...,n - 1. 
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3 Skeletons - classical constructions in billiards 



Before applying the above formalism to construct continuous semiclassical wave functions 
inside billiards B (Fig.l) vanishing on its boundary i.e. satisfying the Dirichlet boundary 
conditions it is first necessary to perform a classical construction consisting of classical tra- 
jectories on which the desired SWF's can be defined. This is just a skeleton i.e. a closed 
set of families of trajectories which forms a base on which SWF's can be constructed. 

The skeleton idea relies on an observation that the short wave packets propagate in 
billiards approximately along the straight lines (rays of the geometrical optics) gradually 
however becoming wider and wider due to unavoidable diffractive effects. Because of the last 
effects even an initially narrow bundle of such rays fills completely a volume admitted by the 
respective boundary conditions. 

Therefore from the one hand the SWF propagates along the classical objects - trajectories, 
but from the other hand it has to be defined on a bundle of such trajectories sufficiently wide 
to gather diffractive effects satisfying nevertheless the rules of the geometrical optics. In 
particular an effect of shadow typical for the geometrical optics should be observed on a way 
of propagation of SWF's. 

Finally a set of all such bundles should be closed and mutually connected to permit us 
performing a construction of the global semiclassical wave function from its pieces defined on 
the separate bundles. 

Below a notion of a skeleton is defined which according to our expectations takes into 
account all the aspects of the semiclassical limit discussed above. 

3.1 Ray bundles and bundle skeletons 

For the needs of this paper we shall assume that the billiards B is classical and according 
to Fig.l it has no holes inside and its boundary dB is a closed curve independent of A and 
given by r = ro(s) = [xq{s), yo{s)] where s is a distance of a boundary point ro(s) measured 
clockwise along dB from some other point A of dB chosen arbitrary, i.e. s{A) = 0. Both 
xo{s) and yo{s) are continues. The curve however consists of a finite number q, q > 1 oi 
smooth arcs or segments Ai, ...,Aq with respective length Li, ■■■,Lq, so that the derivatives 
Xq(,s) and y'^is) are discontinuous in a finite number of points on the segment < s < L 
where L = Li + • • • + is the global length of dB. Both xq(s) and yo{s) are of course 
periodic with the period equal to L. We shall identify the point A with the point beginning 
the arc (segment) Ai. 

Next we define a bundle of rays as a family of trajectories in the following way. 

Let Ak{u, I), Li + • • • + Lfc_i < u < Li + ■ ■ ■ + L^, < I < L^, be an open connected 
piece of the arc (segment) A^ beginning at s = u and having a length I. 

Let further rk{s,t;u,l), u < s < u + I, 0<t, bea family of trajectories given by 
angles 7^(5; u,/), < 7^(5; u,/) < 2tt, at which the trajectories escape from Aig{u,l). The 
angles 7(5; u, I) are smooth functions of s and arc measured with respect to the rc-axis while 
the tangential vectors t(s) = [ ^^^j^-* , ^^^^^^'^ J = [cos /3(s), sin /3(s)] are inclined to the x-axis 
at angles I3{s) (Fig.l). The latter angle can be discontinuous at the points where x'q{s) 
and 1/9(5) are discontinuous. Then the angle ak{s;u,l) = 7fc(,s;?i,/) — /3fc(.s) is made by 
the classical ball momentum p(s;u,l) on the trajectory with the tangent vector t(s), i.e. 
p(s; u, I) • t(s) = pcos akis; u, I). It is assumed that < ak{s; u, I) < it. 

The classical time evolution of the family rfc(s,t;u, Z), u < s < u + l, is therefore the 
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Figure 1: An arbitrary billiards 



following 

rkis,t;u,l) = ro{s) + p{s;u,l)t, ro{s) £ (18) 

where p{s;uj) = [pcos 7^(5; u, /),psin 7^(5; n, /)] satisfies the classical equations of motions 
(11), i.e. p^(s;n, Z) = 2Eq (again we put m = 1 for the billiard ball mass). 

The trajectories (18) define of course the change of variables (x, y) — )• {t, s), in vicinity of 
Ak{u,l),i.e. X = fk{t, s;uj), y = gk{t, s;u, I) with the Jacobean: 

Jkit, s; u, I) = p^'j'kis; u, l)t - p |t(s)| sinafe(s; n, /) = p^j'kis; u, l)t -psinakis; u, I) 

(19) 

since |t(s)| = 1. 

The family of trajectories defined in the above way will be called a bundle of rays 
emerging from the segment Ak{u, I) of A^ and will be denoted by Bi^{u, I) while the trajectories 
themselves will be called rays. 

Note that each bundle by definition is an open set of rays with at most two limiting rays 
as its boundary. If the limiting rays coincide then necessarily the corresponding segment 
ylfc(ti,/) can be closed to the whole billiards boundary. This is the case for example of a 
bundle defined in the circle billiards which rays have the same angular momenta each. 

Suppose the limiting rays of bundles Bk{u,l) and Bki{u\l') coincide on a piece P of 
them. Then by closing the bundles on this common boundary we get a bundle which we call 
compound and denote by i?fc(ti,/) Up B^'^u' ,1'). Such a closing operation will be called 
a composition of the initial bundles. Of course each bundle can be decomposed into two 
others by the reverse operation becoming the composition of the resulting bundles. 

A compound bundle can be composed of course from many bundles and its rays emerge 
then from the sum IJ;. Ak{u, I) of the segments of the composing bundles. If IJ^ Ak{u, I) is a 
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connected piece of the billiards boundary then such a composition and its result a compound 
bundle will be called regular. Other cases of compositions and compound bundles will be 
called singular. 

Since each ray of the bundle Bk{u,l) after some time Tk{s;u,l), ro(s) G Ak{u,l), (dif- 
ferent for different rays) achieves another point of the boundary dB it means that the 
bundle Bk(u,l) maps the segment Ak(u,l) into another piece BAi^{u,l) of the boundary 
dB. In general this map of Ak{u,l) into BAk{u,l) provided by the transformation (18) is 
one-to-one except the caustic points of BAk{u,l) in which JK{Tk{s;u,l),hk{s;u,l);u,l) = 
0, ro(/ifc(s; /)) € BAk{u,l). Here hk{s;u,l), ro(s) G Ak{u,l), realizes explicitly this map. 
If however ^^k{s^u,l) _^ ro(s) G Ak:{u,l), i.e. this map is one-to-one then such a bundle will 
be called regular at the boundary dB or simply regular. 

All bundles considered below will be assumed to be regular. 

By DBk{u,l) will be denoted a domain of the billiards B covered by rays of the bundle 
Bk{u,l) emerging from Ak(u,l) and ending at BAk{u,l). The domain DBk{u,l) is locally a 
Lagrangian manifold on which each loop integral § p ■ dr vanishes. 

While a definition of a bundle is essentially local it can happen that in particular cases 
of bundles or compound bundles they can cover the whole billiards except necessarily some 
pieces of its boundary, i.e. DBk{u,l) = B = D{Bk{u,l) Up Bk'{u',l')). Such bundles or 
compound bundles will be then called global. 

Assume further the boundary dB to be a mirror-like, i.e. reflecting the incoming rays 
according to the reflection principle of the geometrical optics and let y4fc'(u', I') be a piece of 
another arc A^i of dB such that Af^/^u' , I') H BAk{u, I) ^ on which another regular bundle 
of rays Bk'{u', I') = {rfc/(s, t; u', I') : u' < s < u' + I', < t} is defined. 

If on the segment Aki{u',l') fl BAk{u,l) the ray bundle Bk'{u',l') coincides with the 
reflected one we call the ray bundle B^/^u'J') a reflection of the bundle B}.{u,l) on the 
segment mentioned. 

The reflection operation over the bundle Bk{u, I) will be denoted by 11 so that IiBk{u, I) 
denotes the set of all rays arising by the reflection of the bundle Bk{u,l) of dB. 

Consider now a family of the disjoint ray bundles B = (J Bk{u, I), Bk{u, l)nBki{u' , V) = 
ifBk{u,l)y^Bk>{u',l'). 

The family B will be called closed under reflection 11 on the boundary dB if the following 
two conditions are satisfied for each bundle Bk{u,l) G B: 

n 

UBk{u,l) = [j Bj{uj,lj)nUBk{u,l), Bj{uj,lj) G Ai., j = l,...,n 

m 

Bk{u, = U n-Bi(uj, li) n Bk{u, I), Bi{ui, k) G Aj^, i = 1, m 
1=1 

HB = B (20) 

A closed family B' = (J Bk'{u', I') is embedded into a closed family B = [J Bk{u, I) if each 
ray bundle of B' is a subset of some ray bundle of B and each bundle of B contains some 
bundle of B'. 

A closed bundle family is called connected if a unique possibility to represent it by a sum 
of another two disjoint closed bundle families is a decomposition operation done on every 
bundle of the family. 

A closed connected bundle family will be called a Lagrange bundle skeleton or simply 
a skeleton if it cannot be embedded into another closed connected bundle family. 
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Prom now on all considered bundle families will be assumed to be skeletons. 

Let us stress the following four basic properties of skeletons which are of great importance: 

1. each skeleton is complete i.e. none additional bundle can be added to it not destroying 
it connectedness, 

2. each skeleton cannot be decomposed into "smaller" ones (again by the connectedness 

property), 

3. each bundle of a given skeleton is a result of reflections of other bundles of the same 
skeleton (by (20)), 

4. each ray belonging to a skeleton B will never leave B by its time evolution and bounces 
on the billiards boundary. 

Let us now reverse in time all trajectories belonging to B. This operation leads us again 
to some skeleton B"^ which will be called associated with B. 

Bundles of B"^ are obtained simply from the corresponding bundles of B. Namely with 
each bundle Bk{u,l) of B let us associate a bundle B^{u,l) which trajectories satisfy the 
following condition: 

Jk{s;u,l) = Tr + 2P{s)-jk{s;u,l), u<s<u + l (21) 

i.e. these trajectories are just the reflections on dB of the time reversed trajectories defined 
by 7fc(s; u, I) and belonging to Bk{u, I). 

The skeleton B"^ is organized by all bundles B^(u,l) obtained in the above way. 

Needless to say (B^)^ = B. 

Finally let us denote by Db C -B a topological sum of all domains DBk{u,l), i.e. Db = 
DBt(u,i)cB ^^k{u,l)- Define D^a analogously. By the construction of the skeleton B"^ we 
have D^A = Db- 

An useful operation on a skeleton is its reduction which means making all possible com- 
pound bundles of the bundles of the skeleton. Such a form of the skeleton will be called a 
reduced skeleton and denoted by B^. 

Reduced skeletons although useful arc deprived however of many properties of the skele- 
tons. In particular it is typically not possible to construct a skeleton associated with the 
reduced one. On the other hand in the polygon billiards each compound bundle is associ- 
ated with a definite momentum of the billiards ball so that to different momenta of the ball 
correspond different compound bundles. 

If a reduced skeleton B^ contains only global bundles then the skeleton B will be called 
global and if its bundles are all global and regular the skeleton B will be called regular. 
Skeletons which are not regular will be called singular. It then follows that skeletons can be 
global but singular. 

3.2 Skeletons in the phase space. Pseudointegrable billiards 

In the case of the polygon billiards skeletons considered in the phase space are decomposed 
into separate pieces parallel to the billiards plane. Each a piece corresponds to a separate 
compound bundle. Any trajectory of the skeleton visits all pieces (bundles) jumping from 
one piece (bundle) to another when achieving the piece boundary which projected on the 
billiards plane coincides with the corresponding piece of the billiards boundary. 
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One can imagine a continues surface made of the pieces mentioned gluing their respective 
boundaries. Namely, we can get it by gluing points of a piece boundary by which trajectories 
leave the piece to visit the next one with these entry points of the next piece. If a number of 
bundles is finite such a construction leads us to a compact two dimensional surface. If such 
a surface is closed it has then a definite genus. Among the others these are the cases of the 
global skeletons in the rational polygon billiards [8, 13, 23]. 

It will appear in sec. 5 that the skeletons in the investigated billiards can form Lagrange 
surfaces in the phase space which are closed according to the above construction with a 
definite genus but also which are open being a cylinder-like or a Mobius-like bands. 

It will appear also in the next sections that only the skeletons which are regular can 
provide us with GSWF's which are exact solutions of the eigenvalue problems. The singular 
skeletons provides us with singular GSWF's which cannot be exact and some of them show 
typical properties of the supcrscar solutions. 

Billiards for which every possible skeleton has a finite number of bundles are certainly 
distinguished by a possibility of construction of all GSWF's for such a billiards in compact 
finite forms. Therefore it would be reasonable to extend a notion of pscudointcgrability 
introduced by Richens and Berry [8] to such a billiards despite the fact that not all of the 
skeletons corresponding to them can be global so that the continuous Lagrange surfaces which 
construction has been described above can be open, i.e. with boundaries. 

4 SWF's defined on a skeleton 

4.1 BSWF's defined on a bundle 

Consider a skeleton B. On each of its ray bundle -Bjk(«, /) we can now define the following 
pair of BSWF's s; u, I; A), a = ±: 

n{t,s;u,l;X) = 4 ^(t,s;u,0e"^'(^'*+^/>°^'^'=(^'"'')'^^0x-(t,,;^,Z;^ (22) 

where = 2Eq and Xfe(^) ■s; l; A), a = ±, are given by (13) and (15). 

Exactly in the same way we can define a pair '^%j^{t, s; u, I; A), cr = it, of BSWF's on the 
corresponding associated bundle B^{u,l): 

n;k{t,s;u,l;X) = r^^{t,s-u,l)e-'^^^"'-^r^'^'''"'M^ (23) 

It will be also convenient for further considerations to substitute the time variable t by 
the distance variable d = pt and consequently to give the trajectories (18), the Jacobean (19) 
and the wave function (22) the following forms: 

rk{d, s; u, I) = ro(s) + d(s; u, I) 
d(s; u, I) = pt = [dcos7fe(s; u, l),dsm 7^(5; u, I)] 

ro{s) e Ak{u,l) (24) 

and 

Jk{d,s;u,l) = -Jk{t,s;u,l) = — d - smak{s;u,l) 

rois)eAk{u,l) (25) 
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and 

nid, S;U,l;X)= (d, s; U, l)e''Md+i: cos a,(s';u,l)ds') -a ^. I. ^) 

ro{s) e Ak{u,l) (26) 

where Xkid, s; u, I; A) = pxUf^ •^S ^' ^5 ^) (t = ±. 

Nevertheless, for simphcity of notations, the bar over s; n, Z; A) will be dropped in 

our further considerations. 

By the variable d the solutions (15) can be rewritten in the form: 

Xlfi{d,s-,u,l) = xlfi{s;u,l) 
Xk,j+iid, s; u, I) = Xfcj+i(s; u, I) + 



2^ 



^ ( Afc(a,s;u,Z)Xfe,j(a,s;u,Z) + 2 ^ £;j_^+iXfe,„i(a,s;n,0 j da 



m=0 

i = 0,l,2,..., (27) 



where Afc(d,s;'u,Z) = {d, s;u,l)-Ak{d, s;u,l)-J^. ^{d,s;u,l) and Afc(d, s; u, Z) is the Laplacean 
expressed by the variables d and s corresponding to the -bundle. 

^'^((i, s; u, I; A), a = ±, are defined initially in the domain Dk{u, I), D^iu, I) C DBk{u, I), 
of the billiards which boundary dDk{u,l) contains of course Ak{u,l). The remaining part of 
dDk{u, I) is built of the two "limit" rays of Bk{u, I) emerging from the ends of Ak{u, I) and of 
BAk{u, I) if there is no caustic of the bundle Bk{u, I) inside the billiards or by the correspond- 
ing caustic Kk{u,l) = {{fk{Kk{s;u,l),s;u,l), gk{Kk{s;u,l), s;u,l)) : Jk{Kk{s;u,l), s;u,l) = 
0, ro(s) G Akiu,l)}. 

Similarly ^'^.^(d, s; u, I; A), cr = it, is defined in the domain D^{u, I) corresponding to the 
bundle B^{u,l).' 

An important property of the representation (22) is its uniqueness, i.e. for two different 
bundles defined on the segment Ak{u,l) this representation provides us with two different 
pairs of *^(d, s; n, Z; A), a = ±. This conclusion follows from the fact that for A sufficiently 
large the BSWF's are determined only by exponentials and the latter are different at the 
same points {x, y) for different bundles. 

4.2 SWF's vanishing on the bilhards boundary 

Another obvious property of BSWF's ^^(cZ, s; u, Z; A), a = =b, is that they cannot van- 
ish on A}.{u,l) unless XkidjS;u,l]\) vanish there identically. Therefore a wave function 
^'^^'^(x, y; u, /; A) vanishing on Ak{u,l) should be represented in the semiclassical limit by a 
linear combination of at least two BSWF's of the form (22). It is shown in App.A that the 
proper linear combinations have to be the following: 

^l''''{x,y;u,l;X) = ^l{di, sr,u,l; X) + ^^f^id2, S2;u,l; X) = 

Jkhdl,Si;U, /)e"^K*+r cosa,(.';«,0<i.');^-(di, sv, u, I; A) + 

j;.|(d2, sr, u, /)e-^^K<^2-r cosa,(.';«,0d.')^-<^^(d2, s^; u, l; A) (28) 
with the following boundary conditions: 

xms;u,l;X)+x7^ki^,s;u,l;X) = ro(s) G ^fe(«, Z) (29) 
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while the point {x, y) is the cross point of the respective trajectories belonging to different 
bundles, i.e. 

r = [x,y] = rk{di,si;u,l) = ro(si) + di{si;u,l) = 
rA;k{'^-'2i S2; u, I) = ro(s2) + d2(s2; u, I) 
rk{d, s; u, I) G Bk{u, I), rA;k{d, s; u, I) e B^{u, I) (30) 

The vanishing superposition (28) if defined on the bundle By{v',l') suggests that 
'^^"f,,{d,s]u' ,l';X) should be related somehow to the BSWF ^l{d,h^^{s;u,l);u,l);u,l;X) 
defined on the bundle Bk{u, I) which the previous one is a reflection. In the next section this 
relation is established as a condition matching both the solutions. 

4.3 SWF's defined on a bundle skeleton and their continuity 

According to our construction of the skeletons B and B"^ there is a domain Dqb, B D 
Dqb D dB, of the billiards containing the billiards boundary dB in which each point {x, y) 
with its some small vicinities is mapped in the one-to-one way into each bundle of the pairs 
Bk{u,l) and B^{u,l) containing this point. 

Let r = {x, y) <E Dqb be a fixed point of the billiards. Let D{x,y) denote a set of 
all Dk{u,l), Bk{u,l) G B, which contain this point and D^(x,y) is the respective set of 
D^{u,l), B^{u,l) G B"^. According to their definition D};{u,l) G D{x,y) if and only if 
D^iu,l) e DA{x,y). 

SWF's ^^''^{x,y,X) vanishing on the billiards boundary can now be defined on B and 
B"^ in the domain Dqb as follows: 

^^■''"{x, y,X)= J2 in{d{u, I), s{u, 0; u, I; A) + 

Dk{u,l) G D{x,y) 
D^iu,l) e DAix,y) 

^-^,^{d{u,l),s{u,iy,u,l;X)) (31) 

where BSWF's "^'^{d{u,l),s{u,l);u,l;X) and ^^°^((i(u, /), s(n, /); n, A) are defined in Dqb 
and satisfy the condition (29) on Ak{u, I). 

The solutions ^^'~^{x, y, A) and ^'g '~(a;, y, A) coincide if and only if B = B"^. 

The SWF's (31) which satisfy the condition of vanishing on the billiards boundary are 
the most general ones for the skeletons B and B"^ which can be defined in the domain 
Dqb- However the next step in solving the basic problem of energy quantization in the 
semiclassical limit is to make these solutions continuous in Dqb since this property is not 
ensured automatically by (31). ^^''^{x, y, A) are certainly continuous and unique inside each 
bundle contained in Dqb- However if the skeleton B is singular then their bundles have 
their boundaries also inside the billiards area on which '^^''^ {x,y, X) or their derivatives can 
appear to be discontinuous if a point (x, y) crosses these boundaries. They can be also non 
unique if a point {x,y) moves along some closed loops such as the one which is homotopic 
with the billiards boundary, i.e. it can happen that '"^(a;, y. A) do not come back to their 
initial values being continued along such a loop. 

Considering the continuity property of ^'g''^(x,y, A) the following circumstances can ac- 
company in general such bundle boundary crossings: 
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1. two BSWF's ^^(rf, s; u, I; A) and s'; u', I'; A) which enter the sum (31) are defined 
on the bundles Bk{u, I) and Bki{u' , I') which can be composed into a compound bundle; 
and 

2. there is no such BSWF's and the respective neighboring bundles. 

In the first of the above cases the corresponding ^^{d, s;u,l; X), {xk{d, s),yk{d, s)) G 
DBk{u,l), and s'; 1'; A), {xk'{d' , s'),yk'{d' , s')) G DBk'{u',l'), defined by (26) have 

to be identified on the common boundary of the bundles Bit.{u, I) and Byiv! ^ V) together with 
their first derivatives, i.e. 

^^{d,s;u,k\) = ^^,id',s';u',l';\) 
^^Ud,s;u,l;X) = ^^Ud',s'-u',l'-X) 

^^^id,l;u,l;X) = -^H,id',s';u',l';X) 

{xk{d,l),yk{d,l)) = {xy{d\s'),yk'{d\s')) G dDBk{u,l) r^^DBy{u\l') + (32) 

Note that such identifications as the last ones mean that the BSWF's are defined now on 
the reduced skeleton rather then on the original ones. 

In the second case however the corresponding BSWF's ^^{d^ s; n, Z; A), {x{d^ s), y{d, s)) G 
DBk{u, I) or their normal derivatives have to vanish on such a boundary of DBk{u, I), i.e. 

^^id,s;u,l;X)=0 
{xk{d, s), Vkid, s)) G dDBk{u, I) (33) 

or 

^^kid,s;u,l;X)=0 
{x{d,s),y{d,s)) edDBk{u,l) (34) 

that is in such cases ^^{d, s; u, I; A) defined in the bundle B]^{u, I) should satisfy on its bound- 
ary Dirichlet's or Neumann's conditions. 

The last two conditions though necessary seem to look as a little bit arbitrary. However 
we should remember that our calculations are performed in the scmiclassical regime, i.e. 
in the classically allowed regions (bundles) outside which the scmiclassical wave functions 
cannot exist. Physically this means obviously that outside each bundle a corresponding 
piece of the exact wave function represented on the bundle by its respective semiclassical 
approximation has to vanish exponentially (for A sufficiently large) when moving away from 
the bundle. Semiclassically it just means that BSWF's defined inside the bundles have to 
vanish identically outside of them. This condition can course however that the first derivatives 
of the global SWF's (31) can be discontinuous on such bundles boundaries. Just this last 
property differs essentially the semiclassical solutions (31) from the exact ones. If it happens 
we will call such a GSWF singular in contrast to the regular one which is continuous in 
the whole billiards together with its first derivatives. From this discussion it follows also that 
for the latter possibility to happen it is necessary for the skeleton on which *g ''^(x, y. A) is 
defined to be global. 
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4.4 First quantization condition for SWF's 

If ''^(x, y, A) arc made continuous inside Dqb then the uniqueness condition for them 
on each loop lying in Dqb particularly on the ones homotopic with the billiards boundary dB 
if such exists leads us to the first quantization condition which has to be satisfied by these 
two SWF's. 

4.5 Continuing the SWF's ^b''^(x, y, A) over the whole skeletons B and B"^ 
- the global SWF's 

By the formula (31) 'I'b ''^(.x, y, A) are defined in every bundle of the skeletons B and B"^ 
close to the billiards boundary. However by the way of construction of both the skeletons if 
a point (x, y) of the domain Dqb is achieved by a ray of some bundle of B running all the 
time by the domain Dqb then it is also achieved by the same ray inverted in time and being 
a member of a bundle of B"^, i.e. running in the opposite direction. However the second ray 
to achieve the considered point {x,y) G Dqb has in general first to leave the domain Dqb 
crossing its boundary in several points in order to come back to it. 

The same note is valid for rays contained in the skeleton B"^. 

Therefore {x,y,X) defined by the formula (31) can be continued from a point (x,y) 
of the domain Dqb into another such point of Dqb along rays contained in the skeletons B 
or B"^. If Dqb cannot be equal B, then such a continuation meet as necessary caustic points 
which have to be avoided somehow. If we do that however we will achieve again points of the 
domain Dqb and naturally the continued solutions and the solutions defined by (31) have 
to coincide. This coincidence formulate the second quantization condition which both the 
SWF's *B''^(x,y, A) have to satisfy. Such a coincidence is achieved by identifying each term 
of the sum (31) with the corresponding term of the continued ''^(x, ?/, A). Anticipating the 
results of App.B of [22] the corresponding identification should be done as follows. 

1. Let Bt;{u,l) be a reflection of the bundles Bk^{u\,l\),Bk^{u2,l2), ■ ■ ■ : Bk^{un,ln) satis- 
fying (20). Let ^l'^°'^\d, s; ui,li; A), (d, s; U2, ^2; A), . . . , 'i>'^^'^\d, s; u,,, In, A) de- 
note the SWF's defined in Dqb and continued on the respective bundles -B/jj (tti, ^i), 
Bk2{u2,l2), ■ ■ ■ ,Bk„iun,ln) again to Dqb- Let further ^'^|^(d, s; u, Z; A) be defined in 
Dqb on the bundle B^{u, I) while *^(d, s; u, I; A) in Dqb on the bundle Bk{u, I) being 
both related by the boundary condition (29). 

Then we make the following identification of BSWF's: 

'^A^k((^,Hs;uj,lj);u,l;X) = '^l'^°'^\D{s]Uj,lj) - d,s;uj,lj; X) 

ro{h{s;Uj,lj)) G Ak{u,l) Ci BL{uj,lj) 
ro(s) e L{uj,lj) 

ro{h{s;uj,lj)) = ro{s) + 'D{s;uj,lj), j = l,...n (35) 

where D{s; u, I) = |D(s; u, l)\ denotes the distance between the points ro(/ifc(s; u, I)) and 
ro(s) of the boundary dB. 

Similarly 

2. Let B^{u,l) be a reflection of the bundles Bf^{u[,l[), Bf^{u2,l2), ■ ■ ■ , Bf^{u'j^,l'^) sat- 
isfying (20). Let ^'2^'\d, s; u[,l[; A), ^T^^, s; 4, A), . . . , ^Z'^^, s; u'^, 4; A) 
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denote the SWF's defined in Dqb and continued on the respective bundles Bj^{u\,l'i), 
B^{u2,l2), ■ ■ ■ , Bjl^(u'^,l'^) again to Dqb- Let furtlier ^^'^(d, s; ti, A) be defined in 
Dqb on the bundle Bk(u, I) while *^.j;.(ci, s; u, I; A) in Dqb on the bundle B^(u, I) being 
both related by the boundary condition (29). 

Then we make the following identification: 

'^^''{d,h{s;ui,li);u,l;X) = *^'.^°"*(_D(s; Uj, /j) - d, s;ui,li; X) 

ro{h{s; Uiji)) e Akiu, I) n BL{ui, k) 
tq{s) G L{ui,li) 

ro{h{s-Ui,li)) = ro{s) + 'D{s;ui,li), i = l,...m (36) 

3. Meeting the caustic points the BSWF's ^'^(d, s; u, I; A) and "^^/^id, s; u, I; A) avoid them 
by fixing s and moving on the complex c/-plane from above the points for cr = + and 
from below them for a = —. 

The conditions (35) - (36) allow us to define *b ''^(z, y, A) as given by (31) in every 
point of the domain Db, i.e. in the domain classically allowed when moving on the skele- 
ton B and to rewrite the sum in (31) representing ''^(z, y, A) globally by the terms of 
^'j^{d{u,l), s{u,l);u,l; X) or by the terms of "^^"'i^{d{u,l), s{u,l);u,l; X). Namely, the global 
SWF's (GSWF) '"(a^' y> ^) are given by: 

*B 2/> A) = E nid{u, l),s{u, iy,u, I; A) = 

DB^{u,l)&D{x,y) 

^^,{d{^,l),s{u,iy,u,l;X) (37) 

where D{x,y) and DA{x,y) denote now the respective sets of DBk{u,l) and DB^{u,l) con- 
taining the point (x,y). 

The sums in (37) contain all BSWF's '^^.{d, s;u,l; X) and *^^^((i, s; -u, Z; A) which can be 
continued to this point by the corresponding domains DBk{u,l) and DB^{u,l). 

Let us note that the forms (37) of the GSWF's allow us in fact to define them on the 
reduced form of the skeleton B rather then on the skeleton itself. This possibility permits 
to reduce substantially number of terms in sums (37) . 

Rewritten in terms of the x-coefficients Eq.(35) gives: 



XA%(.d,h{s;uj,lj);u,l;X) 



aiXpSk{uj,lj) 



dh{s; Uj, Ij) 



ds 



Xk.°^\D{s- Uj,lj) -d,s; Uj,lj; A) 



rs rh(s;uj,lj) ^ ^ 

6k{uj,lj) = D{s;uj,lj) + / COS {s ;uj, I j)ds — / cos ak{s ■,u,l)ds 

Juj Ju 

ro{h{s; Uj , Ij)) G Ak{u,l) Ci BL{uj,lj) 
ro(s) G L{uj,lj) 
ro{h{s;uj,lj)) = ro{s) + D{s;uj,lj), j = l,...n (38) 

Note that 6k{uj,lj) in the above formula is s-indcpcndcnt (see App.B). Due to that and 
due to the properties (16) and (17) the rhs of (38) satisfies (12) as it should. 
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Putting d = in (38) and taking into account (29) we get: 

Xk{^,h{s;uj,ljy,u,l;X) = -x^^kiO, h{s;uj, I j);u, I; X) 



-Vae 



dh{s] Uj, Ij) 



ds 



^a^cont ^D(^s-Uj,lj),S-Uj,lj-\) = 



-r/,e«(«-'^)x^f "*(0, h{s- Uj,lj)-u, /; A) 
ro(/i(s; Uj,lj)) G Ak{u, I) n BL{uj,lj) 
ro{s) G L{uj,lj) 

rQ{h{s;Uj,lj)) = ro(s) + 'D{s;Uj,lj), j = l,...n 



(39) 



The BSWF's *^(d, s; n, Z; A) defined on bundles of B and *^(d, s; -u, /; A) defined on re- 
spective bundles of B"^ are related with each other by the boundary conditions (29) and by 
matching conditions (36)- (39). 

It is clear that the conditions (39) have to determine also the factors x^(s;?x, Z;A) 
for all the bundles Bk{u,l) which are the "initial" conditions for both x^(d, s; n, Z; A) and 
X%kid, s; u, I; X) in the recurrent formula (15), i.e. x^(s;u, Z;A) = Xfe (Oi s; u, Z; A) = 
— X^.fc(0, s; u, /; A). Nevertheless these conditions cannot be given arbitrarily. Just opposite 
all Xkis', u, I] A) have to satisfy (39) in a selfconsistent way. 

The formulae (38) and (39) define the conditions which the SWF's Xfc(c^) H^'i "^jJj)'-: ^5 ^) 
should satisfy when bouncing from the billiards boundary. Nevertheless this condition can be 
specified additionally with respect to its factors. Namely, taking their large A-limit we get: 



XkAKs\Uj,lj);u,l) 



dh{s; Uj, Ij 



ds 



Xk,r+iih{s;uj,lj);u,l) = -r]^ 

rD{s;Uj 



dh{s; Uj,lj] 



ds 



Xki,r+l{s;Uj,lj) + 



2pJo 



5x^.^^(a,s;«j,Zj)) + 



2^£;r_j+iXfe^,;(a,s;Uj-,Zj) da\ 
1=0 ) ) 

r = 0,1,2,..., 



(40) 



The above equations should be satisfied on each bundle Bk{u, I) of the skeleton B. 

The first of the equations (40) should determine the classical quantities, namely the 
skeleton B and the "classical" energy Eq = ^p^ and by them define the JWKB approximation 
of the SWF's. Namely: 



*r^^-(x,y,A) = 



,JWKB:a 



J2 J Hd{u,l),s{u,l))e 

DB{u,l)GD{x,y) 



DB{u,l)eD{x,y) 
aiXpl d{u,l)+ /j'"''^ cos a{s';u,l)ds' 



''^{d{u,l),s{u,l);u,l) = 



xUdiu,l),siu,iy,u,l) (41) 



The remaining equations determine quantum corrections to the "classical" ones involved 
in (41). 
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However it is easy to note that for the selfconsistency of the equations (40) it is necessary 
for the exponent e'^i^P^ki^j ]-)g independent of A, i.e. we have to have on each bundle 
Bkiu,l) of B: 

Xp6k{u,l) = (pkiuj) 

Bk{u, I) CB (42) 

where 5k{u,l) is given by (38) and (f)k{u,l) is a A-independent constant. 
The equations (42) have to define the energy Eq = ^p^. 

Taking into account the last conclusions we get the following final set of the recurrent 
quantization conditions: 



Xp6k{u,l) = (pkiuj) 



ai<j>kAuj,lj) 



dh{s; Uj,lj) 



Xk,r+\{Ks\Uj,lj)\u,l) = -rjae-^ 

(71 f^i^''^j>h) 

2^ 



ds 

dh{s; Uj,lj 



ds 



X%fl{s;uj,lj) 



r = 0,1,2,..., (43) 



together with: 



Xkfl{d,s;u,l) = Xk,o{s;u,l) 

Xk,r+lid, S; U, I) = Xk,r+l{s\ U, I) + 



r = 0,1,2,..., (44) 

Let us note finally that if B 7^ B"^ then energy levels corresponding to the skeleton B have 
to be degenerate. This conclusion follows easily from the form of the quantization conditions 
(43)-(44) and (32)-(33) showing that the complex conjugations of *g ''^(x, y. A) satisfy also 
these conditions with the same semiclassical energy E. The two corresponding solutions are 
of course ^'g '^(x, y. A). 

4.6 Finite and infinite bundle structures of skeletons. The last quantiza- 
tion condition 

The constructions of skeletons and SWF's in billiards performed in sec. 3-4 describe com- 
pletely the energy quantization problem in the semiclassical approximation. 

For a given billiards however there can be skeletons with a finite number of bundles 
as well as with an infinite one. The semiclassical quantization procedure described in the 
previous sections seems to be easily applied to the finite bundle number skeletons. Namely 
in such a case following a trajectory starting from a bundle Bk{u, I) we have to approach the 
same bundle after a finite number of bounces. The corresponding semiclassical wave function 
propagated by the skeleton has therefore to come back to its initial form achieving again the 
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initial bundle. This condition closes essentially the process of quantization formulated in the 
previous sections. The respective conditions are of course the following: 

exp I Yl ^kiu,l) I n(-?7aj = 1 

\Bfc(it,OeB / k=l 

Xk'^\D{s,s),s;u,l;\) =Xk{s\u,l;\), s,s'eBk{u,l) (45) 

where n is a number of bounces and D{s,s') is the global distance passed by the billiards 
ball along the investigated trajectory. 

Skeletons with a finite number of bundles are typical for the billiards with the integrable or 
pseudointegrable motions. Nevertheless they can be found also as particular cases of motions 
in chaotic billiards as well. 

The cases of skeletons with an infinite number of bundles arc clearly much more difficult 
for investigations. Such skeletons should be typical for chaotic billiards. 

According to its definition a bundles Bk{u,l) can bifurcate after the reflection by the 
billiards boundary into many diff'crcnt subbundles, i.e. parts of other bundles having their 
beginnings also partly on the arc BAk{u, I). In fact a general behaviour of a skeleton in such 
chaotic cases should not differ essentially by its chaotic complexity from a chaotic trajectory 
reminding however rather a gigantic road-knot with infinitely many viaducts spanning the 
billiards boundary on which the billiards ball moves. It is obvious that if they exist their 
identification seems to be not an easy task. 

Nevertheless the rule (46) can appear to be useful also even in such cases. This is because 
a ray beginning with a bundle Bk{u, I) can come back to it even arbitrarily close to its initial 
starting point on Ai.{u,l) (according to the Poincare theorem) not repeating its way. But 
this is enough for writing the " last quantization condition" (46) where the sum goes now over 
all bundles of the skeleton passed by the ray. 

In the next two sections we shall focus on the finite number cases of bundles in skeletons, 
i.e. applying this procedure to the simplest well known cases of the polygon billiards not 
avoiding however billiards with chaotic motions such as the Bunimovich one. 



5 The rational polygon billiards 

A two dimensional rational polygon billiards are distinguished by their pseudointegra- 
bility [8]. As we have discussed it in sec. 3 a phase space corresponding to a motion in such 
a billiards on a given skeleton consists of a finite number of pieces parallel to the billiards 
plane and orthogonal to the two momentum axes and corresponding each to the compound 
bundles the reduced skeleton. In the case when the corresponding skeleton is regular then 
by the gluing procedure described in sec.3.2 one can get [8, 23] a two dimensional compact 
closed surface with a genus g given by: 

where N is the number of the compound bundles, n is the number of the polygon vertices, 
and TT^ with integers pk, Qk relatively prime is the angle enclosed A;-th vertex, k = 1, ...,n. 

In other cases of the skeletons developed in the rational polygons one gets surfaces which 
do not provide us with closed surfaces in the phase space, i.e. such skeletons are singular. In 
particular such singular skeletons are developed by periodic trajectories. 
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Note that it is the polygon skeleton property that if it contains at least one periodic 
trajectory then all trajectories of such a skeleton are also periodic. In such a polygon periodic 
skeleton there are always two (and no more) periodic trajectories each of which starts from 
some vertex of the polygon and runs to another one. These two periodic trajectories have 
been called by Bogomolny and Schmit [17] as singular diagonals (SD's) while the skeleton 
itself as the periodic orbit channel (POC). Therefore each periodic skeleton is defined by two 
SD's. 

A convenient way of representing motions in a polygon billiards can be obtain by unfolding 
the polygon by its repeating reflections in its sides on which the trajectory reflections are 
performed. A frequently complicated pattern of the real trajectories takes then a simple form 
of parallel straight lines on such unfolded polygons. 

While a triangle is the simplest polygon its billiards properties are in general not as such. 
A motion in rational triangles can be intcgrable if g = 1 for them so that for the triangle 
angles tt^, i = 1,2, 3, as it follows from (46) we have to have: 

Pi = l, ?: = 1,2,3 
111 

- + - + - = 1 (47) 

qi 92 93 

Several obvious solutions to (47) give the following triangle angles (^, |, ^), (|, j, ^) and 
(f ' f ' f ) integrable cases. 

In fact it is rather a rectangular billiards and its variations which we call broken rectangle 
billiards which seems to demonstrate sometimes in a spectacular way most advantages of 
the skeleton approach developed in sec. 3. Therefore we will firstly consider the cases just 
mentioned. The cases of the equilateral triangle and the pentagon billiards will be considered 
next. 

5.1 The rectangular billiards 

Consider therefore the rectangular billiards shown in Fig. 2. This billiards is the canonical 
example of the energy quantization problem because of its easiness to be solved by the variable 
separation method. According to Fig.2 the well known solution to the problem is given by 
the following two equations: 

XpxCi = mvr 
Xpyb = mr 

m,n= 1,2,3,... (48) 



giving the energy: 



„ 1 r, 1 r, vr^ / \ , , , 

^=?' + ??=2A^U+(? (49) 



and being the result of the following form of the (non-normalized) energy eigenfunctions: 

*(x,y) = 4:sm{\pxx) sm{Xpyy) = 

^XpxX-Xpyy _|_ ^-XpxX+Xpyy _ ^XpxX+Xpyy _ ^-XpxX-Xpyy ^qq^ 

which have to vanish on the billiards boundary. 
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• B(0,a;a) 

B(0,a;n-a) 

■ B(a,b;7t-a) 

■ B(a,b;7i+a) 



.B(a+b,a;ji+a) 

- B(a+b,a;27i-a) 

- B(2a,b;2re-a) 

- B{2a,b;a) 



Figure 2: The eight ray bundles of the generic regular skeletons B and B (= B) in the 
rectangular billiards and the one {Bg U Bi) of the four compound bundles. The latter form 
a torus in the phase space. 



Of course one can always put px = p cos a, py = psina where a, < a < ^t^, is the angle 
by which the momentum p is inclined to the x-axis when the billiards ball reflects from the 
side Ai. Therefore the classical trajectory angles of the billiards ball are quantized according 
to the formula: 

n a 

tana = — — 
m 

m, n = 1, 2, 3, ... (51) 

Let us note that the cases a = 0, ^vr are excluded by the solutions (50). 

Let us note further that the set S of all pairs {m,n), m^n = 1,2,3,..., defining the 
eigenfunctions '^rn,n{x,y) can be divided into disjoint subsets Smo,no each of which contains 
a pair (mo,?T-o) where mg and no are relatively prime and all its multiples {kmQ,kno), k = 
1, 2, 3, and nothing more so that S = IJ^^j Ilmo,no- It is clear that all the points of Smo,no 
lie on the straight line y = tanao^; with tanao = ! i-e. all the states ^kmo,knoix,y), k = 
1, 2, 3, are related in the rectangular billiards with a family of classical trajectories which 
are inclined to the x-axis by the angle qq. 

5.2 The rectangular billiards skeletons built by non-periodic trajectories 

To perform semiclassical calculations corresponding to "generic" skeletons let us consider 
a skeleton shown in Fig. 2 containing, by assumption, only nonperiodic trajectories. According 
to the description of the previous section there are four "smooth arcs" in the rectangular 
billiards, i.e. the four sides of the rectangle Ai,...,A^. Since the absolute values of the 
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momentum components Px,Py are the integrals of the classical motion inside the billiards 
respecting elastic law of bouncing then all bundles which should be taken into account are 
defined by a single angle a, < a < ^tt, which are made by the rays of the bundle Bi = 
Bi{0,a;a) with the x-axis. 

Choosing the case of the angle a shown in Fig. 2 the remaining seven bundles of the 
skeleton B shown in this figure are: B2 = -Bi(0, a;7r — a), -83 = B2{a,b;7r — a), B4 = 
B2{a, 6; TT + a), B^ = B^{a + 6, a; tt + a), Bq = B^{a + b, a; 2'k — a), B^ = ^4(20 + 6, b; 27r — 
a), Bg = B4(2a + b,b; a), i.e. the parameter s introduced in sec. 3 is counted anticlockwise 
starting from the point (0, 0) of Fig. 2 (and having negative value if measured clockwise). The 
bundles B2k-i, B2k, k = 1, ...,4, are defined on the respective sides Aj^, k = 1, ...,4, of the 
billiards, i.e. on Ai = Ai{0,a), ^2 = ^2(0,6), A3 = As{a + b,a), A4 = ^4(20 + 6, 6). 

The skeleton coincides exactly with B in the case of generic skeletons in the rectangular 
billiards, i.e. the corresponding energy levels cannot be degenerate. 

Let us note that a number of bundles in the skeletons is obviously independent of a choice 
of a, i.e. it is always equal to eight if < a < ^tt. 

Thecorresponding Jacobean factors of A), q = 1,...,8, are J~2 (^d, s) = (— sinafe)~2, 

Ok = a, s £ Ak, k = 1,3 and = ^vr — a, s € A^, k = 2,4, i.e. the Jacobeans are constant 
but discontinues. Therefore they will be included into the x-factors contained in the SWF's. 

We can now make use of the fact that from the sixteen BSWF's ^^{d, s. A), q = 1, 8, 
we can first select only eight of them with the positive signature since the negative signature 
solutions have to coincide with the respective positive signature ones. Next since each pair 
s. A), 'i>2q^i{d, s, \) q = 1,...,4, of these solutions has to coincide on the common 
boundary of the respective bundles B2q and B2q+i, q = 1,...,4, then we can define the 
solutions on the respective compound bundles to get in this way only four BSWF's, namely: 

^^{d,s,X) -b<s<0 

^f{d,s,X) 0<s<a 

^2{d,s,X) 0<s<a 

^f{d,s,\) a<s<a + b 

^'4(^,5, A) a < s < a + b 
^f{d,s,X) -a-b < s < -b 

^'g(<i, s,A) -a-b < s < -b 
^'|(d,s,A) -6<s<0 

and also the respective four compound bundles: 



Bi 


= -Bs,! 


= 58 U Bi 


B2 


= -62,3 


= -B2 U S3 


B3 


= -64,5 


= 54 U 


B4 


= -66,7 


= Be U 



(53) 

on which the four solutions (52) are defined. 

Therefore the reduced skeleton B^ contains four compound bundles (53). Note that 
each compound bundle Bk, A; = 1, 4, is regular so that the "generic" rectangular billiards 
skeleton B is also regular. 



*+(d,s,A) = 

^+{d,s,X) = 

^tid,s,X) = 

^tid,s,X) = 
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X a 



Figure 3: The four rays of the four compound bundles and the corresponding four solutions 
'^f{di,Si), i = 1,...,4, meeting at the point {x,y) 



The above compound bundles are shown on Fig. 3 where in each billiard point four their 
rays are met and the four solutions (52) are superposed to get GSWF, i.e. 

4 

M/'^^(x,y) = ^§+(4,Sfc,A) (54) 
fc=i 

Since in our further considerations we will work exceptionally with the solutions (52) we 
will drop the signature of these solutions as well as the tilde mark for a convenience. Then 
assume the following standard forms for ^k{d, s,X), k = 1, 4, : 

^k{d, s, A) = e^WApsco.a,^^(^^ ^) (55) 

where is the angle the momentum of the ray of the compound bundle Bk makes with the 
corresponding side of the rectangle measured anticlockwise. 

The solution ^k{d, s,X) is defined on the compound bundle Bj. which rays start from the 
sides and so that the variable s is measured from the left end of the corresponding 

side Ak, k = 1,...,4. For a given ^^((i, s,A) s is then positive on Ak and negative on Ak-i 
where ^k{d, s, A) is also defined. Since the constant Jacobean factors have been included into 
X-coefficients the coefficient Xk{d, s, A) is continuous on the sides A^-i U Ak, k = 1, 4. 

Let the colours of rays corresponding to the particular compound bundles B^ denote also 
colours of these bundles. Then unfolding the skeleton of Fig. 3 onto the plane a motion of the 
billiards ball which begins with the rays of the bundle Bi is limited by the stripe bounded 
by the two parallel (thick black) lines shown in Fig. 4 and are performed along the straight 
line. It is seen clearly on the figure that this motion is just the scattering of the skeleton 
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bundles on the (white) vertices of the rectangle so that each bundle is scattered into the two 
neighbour ones with the exception of its single ray which crosses the vertex. The latter ray is 
scattered back into the third remaining bundle. No one of the straight line rays crosses the 
rectangle boundary at the same point and the crossing points of each ray are dense on each 
rectangle side. The ray of Fig. 4 which starts at the point with the coordinate s on the figure 
is also shown in folded way on Fig. 5. 

The GSWF (54) has to vanish on each side of the rectangle and this condition exhausts 
all the conditions it has to satisfy. However, a particular form of the corresponding conditions 
depends on a choice of points on the rectangle billiards boundary even for the same bundle. 
Therefore let us choose for writing these conditions the four first points of the ray (including 
the starting point) shown in Fig.5. as convenient for our further considerations. As it follows 
from the form (54) of the solution and from Fig.4 and Fig.5 the corresponding conditions are: 

g-apssina s, A) + e^^P(«^°«"-^«i°")x2 ( -a - s cot a, x] + 

V sm a / 

b + s , ,\ 

, a — [0 + s) cot a, A + 

sma / 

^iXp{b+s)sina^^^Q^ 6 + s, A) = (56) 



,iXp{a cos a+(&+s) sin a) 




^iXpi — hs cot a cos a) f b -{- S 



Xi -. ,5, A + 

sma 



gap(6sina+(a-(6+5)cota)cosa) / _^ + s) COt a - a, A^ + 

Vsma / 

gap((a-(6+s) cota) cosa)^g(-Q^ ^ - + S) COt a, A) + 

gUp(6+s)cotacosa^^(0, -(6 + s) COt a, A) = (57) 



gap(26+s)cotacosa^^(Q^ (26 + s) COt «, A) + 
„iAp(acos a-(2fe+s) cot a cos a)^^(Q^ (26 + s) COt a - a, A) + 



ap(acosa+6sina-(26+.)cotacosa)^^ ( a + ^) '^"t « ^ 35 + ^ _ ^ tan a, A^ + 

\ COS a / 

e^^f(iiliJ+(''+^)^°t"'=°«")x4 f^,-(6 + s)cota,A^ =0 (58) 

Vsma / 

gap(^-(2.+.)sin.)^^ f^^i?^±^l^,(26 + .)cota,A') + 



cos a 

^i\p(a tan a sin a— (26+s) sin a) ^ 



X^iO, a tan a — 26 — s, A) + 

gap(-atanasina+(36+s)sina)^^(0^3^ + S - Otana, A) + 



gap((36+.) sin a-a tan a sin a+a cos a) / 36 + g atana ^ ^ a - 2a, A^ = (59) 

\ sin a / 

where according to our convention s is negative being measured from the left end of the side 
Ai of the rectangle. 
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Figure 4: The unfolded motion in the rectangular billiards 




Figure 5: The solution ^i{d,s) being carried by the successive bundles Si — )• B4 — )• — )• 
i?2 — )• -B3 — -B4 — )• -Bi — ••• • The corresponding ray is shown also unfolded in Fig. 4 
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The last equations reduce to the following ones: 

XI (0, s, A) + e^^f«'^°^"x2 -a - s cot Q, a) = 

V sin a / 

gUpacosa^g ft^^a -{b + s) cot a, a) + X4(0, 6 + s, A) = (60) 



^iXpbsina^^ I t^^s, A ) + X4(0, -(6 + s) COt «, A) = 



sm a 



gUpftsina^^ (26 + s) cot a - a, a) + X3(0, a -{b + s) cot a, A) = (61) 

Xi(0, (26 + s) cot a, A) + e'^P''«'°"x4 f — -(^ + s) cot a, a) = 



V sin a ' 



X2(0, (26 + s) cot a - a, A) + e^^^^^^^^xs ( - — (^^ + ^)^°^"' , 36 + s - atana, A^ = (62) 

V cos a J 



gUpacosa^^ ( a — (26 + s)cota,A^ +X2(0,atana-26-s,A) = 
V COS a J 

X3 (0, 36 + g - a tan a, A) + e^^^" cos + g - Q tan a ^ ^ ^ a - 2a, a) =0 (63) 

V sm a J 

The first two equations of (61) and (62) can be "solved" in the first order in A with the 
help of (15) to get: 

Xi,o((26 + s) cot a) = e^'^P'^'^^xifii'^) (64) 

The point s = (26 + s) cot a lies now on the side Ai of the rectangle. It is achieved by the 
ray leaving the starting point s of the side ^4 of the rectangle and passing the distance . 
Since the coefficient Xi,o(s) does not change propagating along the ray we have to have: 

Xi,o((26 + s) cot a) = Xi,o(s) (65) 

so that 
and 

Ap6sina = n7r, n = l,2, .... (67) 
Similarly, choosing another propagation path we get: 

g2iApacosa _ ^ ^gg-j 
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and 



Xpa cos a = miT , m=l,2, .... (69) 

On the other hand the equation (65) is of great importance because it shows that Xi,o(s) 
is defined in different points of the sides Ai , ^4 by its value estabhshed in some definite point 
of these sides and by the propagation procedure defined by (14) and (15). Since however a 
propagated ray reflects consecutively on the boundary in points densely distributed on it this 
initial value is also propagated densely on the boundary. Therefore to get the coefficient as 
a continuous function of s we have to put it a constant. Let it be equal to one. 

The next order term propagates according to the formula: 

Xi,i{d,s) = xi,iis) + '-yd (70) 

where Xi,i(*) is the initial value of the term on the boundary. Therefore, since distances di 
measured along the propagating ray of the consecutive boundary points Sj by which the ray 
is refiected are distributed on the boundary irregularly but densely values of Xi,i(0)'5j) = 
Xi,ii^i) ~ Xi,iidi,s) in these points have to change discontinuously. Therefore to maintain 
the continuity property of Xi,i(^) o^i the boundary we have to put £^1 = in (70) so that 
Xi,i{d,s) is again constant on the boundary and in consequence independent also of d. 

Quite similarly we can argue that also the remaining terms have to be constant as well 
as all the other terms of the energy semiclassical series have to vanish. 

In this way we get finally from (61)-(64): 

Xl,mnid,S,X) = 1 
X2,mn{d,S,\)^{-ir^' 
X3,mn{d,S,X)^{-ir+^ 

X4{d,s,X)^i-ir+' (71) 

and 

1 



Emn = Eo,mn = ^Pmn = ^ (^"^ "P" J ' 1'n,n = 1,2,3, ... (72) 

Choosing therefore the point {x,y) of Fig.3 for the SWF (54) we get: 

4 4 

^mni^, y) = Yl ^lidhsi) = J2 e'^P'-'^'^^+'^P'-'^f'^''^+'l" (73) 
1=1 1=1 

where di and fi{si), I = 1, ...,4 should be calculated from the relations (see Fig.5): 

X = xo{si) + di cos a{si) 
y = yo{si) + di sin a{si) 

and the phases are defined by (71). 
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However making use of the independence of the phase integral /^q Pxdx + Pydy of the 
integration paths we get for the particular terms in the sum in (73): 

giApd2+iAp/2(s2)+j(n+l)vr _ ^i\(px{a-x)+pyy+(n+l)-K) _ 
^iXpd3+i\pf:i{s:i)+i{m+n)'K _ ^i\(px(a-x)+py{b~y)+(m+n)iT) _ 
giApd4+iAp/4(s4)+i(m+l) _ ^i\{jpxX+py{b-y)+(m+l)-K) _ 

where = Pmn cos amn, P=Pmn sin amn- 
Finally: 

(x, y) = e^^^'^^+Pyy) — ^^K-P^cx+Pyy) _|_ ^iK-P^cX-pyy) _ ^iX{pxX-pyy) _ 

-4sin(pj;a;) sm{pyy), m, n = 1, 2, 3, ... (76) 

reproducing in this way the exact result (50). 

Let us note that the set of all SWF's (76) is again complete. 

Of course it is not surprising that the scmiclassical calculations performed above repro- 
duce the exact result (50) since it follows from its form that it represents simultaneously its 
scmiclassical expansion. However to get the result (76) we have had to assume that skeletons 
considered had to be "generic", i.e. they were constructed of non-periodic trajectories. Just 
this assumption allows us to use the arguments of dens distributions of values of the SWF to 
establish its value on the boundary. In fact no other argument exists to get such a conclusion. 
On the other hand such an argument cannot be invoked in the cases of skeletons which are 
built of periodic trajectories so that such cases of skeletons must be considered separately. 

Finally let us conclude that: 

1. each non-periodic skeleton in the rectangular billiards is regular and equivalent in the 
phase space to a two dimensional torus; 

2. the GSWF (76) is obviously regular and provides us with the exact solution to the SE; 

3. the GSWF (76) as well as the corresponding energy levels coincide with their JWKB 
approximations since the corresponding semiclassical series for the SWF's and energy 
levels abbreviate on the zeroth term; and 

4. it is the dense distribution of the skeleton rays in the configuration and the phase spaces 
which causes the semiclassical series abbreviation mentioned. 

5.3 Skeletons built by periodic trajectories 

Periodic skeletons in the rectangular billiards can be easily realized since bundles of such 
trajectories with a given period are defined by one of its members which starting from a 
vertex of the rectangle has to be reflected in other vertices to " finish" its motion in the initial 
vertex. Such a leading particular periodic trajectory will be called a singular diagonal (SD) 
after Bogomolny and Schmit [17]. All other trajectories which make the same angles with 
the corresponding sides of the rectangle as this particular SD does are then also periodic. 

In the rational billiards periodic skeletons are defined always by two such SD's. In the 
rectangular billiard these two SD's are symmetric with respect to each other in a sense that 



_ f,iMPxX+pyy) 
_QiH-PxX+pyy) 
, ^iH-PxX-pyy) 
_JMPxX-pyy) 



(75) 
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Figure 6: Three unfolded skeletons defined by the corresponding pairs of SD's. The skeletons 
// and /// are singular. Every skeleton is the cylinder-like Lagrange surface in the phase 
space. 



if one of them crosses some two different vertices of the rectangle the second one has to cross 
the remaining two. 

On Fig. 6 and Fig. 7 there are shown five cases of skeletons defined by pairs of such SD's 
represented in their unfolded form (the left picture) and in their real form in the billiards 
(the right pictures). Each pair of SD's defining each skeleton are visible in the unfolded form 
of the skeletons as two parallel straight lines. Each skeleton is a stripe bounded by such two 
SD's. A general property of each such a stripe is that all the rectangle vertices related with 
the stripe lie on its boundary, i.e. on its two SD's. 

Single periodic trajectories are shown also in each skeleton case in the figure being parallel 
to the SD's defining skeletons. In the billiards (the right pictures) these periodic trajectories 
are of course closed. The skeletons on the figure have forms which are typical, i.e. infinitely 
many others differ from these on the figure by a number of refiections of SD's on the rectangle 
sizes. 

5.3.1 Bouncing ball skeleton and the corresponding regular GSWF's 

We will construct GSWF's on these skeletons with the same rules as formulated earlier. 
Let us begin with the case of the skeleton numbered by I in Fig. 6 and shown in Fig. 8. 
This bouncing ball skeleton contains only two bundles Bi and -63 - the first one with its rays 
directed up and starting from the side Ai and the second -B3 with rays directed down starting 
from the side ^3. The skeletons B"^ and B^ are identical with B which is of course regular. 

For these particular cases of bundles rays for both the bundles will be positioned by the 
same parameter s measuring a distance of a ray from the y-axis along the corresponding sides 
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Figure 7: Another two unfolded singular skeletons defined by the corresponding pairs of SD's 
with the same cylinder-like Lagrange surfaces in the phase space. 
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Ai and A^. Therefore for the corresponding BSWF's: 

*3(6 - d, s, A) = e'^P^'-'^'^xsib - d, s, A) 

0<d<6, 0<s<a (77) 

For the coefficients Xk{d,s,X), A; = 1, 3 it is assumed as usually that they propagate along 

the rays of the bundles continuously and this their property is not influenced by reflections 
of the rays on the boundaries. Therefore we have to accept also that they are periodic with 
respect to the d-variable with the period equal to 2b. 
For the GSWF y, A) we have therefore: 

^'''{x,y,X) = ^i{y,x,X) + ^3ib-y,x,X) (78) 

together with the following boundary conditions on the sides Ai and A3 respectively: 

Xi{0,x,X)+e'^P'>Xsib,x,X) = 

e'^P'xi{b,x,X)+X3{(i,x,X) = (} (79) 

so that: 

Xi,o{^) = e^'^^'xi,o{^) (80) 

As previously we conclude that: 

Xpb = mr, n= 1,2,3,... (81) 

Now the corresponding boundary conditions for ^k{y,x, X), k = 1,3 on the sides A2 and 
A4 give: 

Xfc(y,o, A) = xk{y,a,x) = o 

A; = 1,3 (82) 

Therefore in the zeroth order we have: 

Xi,o(0) = Xi,o(a) = (83) 

Next let us invoke the second of the equations (15) and the periodicity of xi(y,x, A) to 
get in the considered case for the second order term: 

Xi,i(26,x) = xi,i{0,x) = xi,i{0,x) + ^ (^^^^ + ^E,xi,o{x)^ (84) 

so that 

f^l^ + 2E,xi,o{x) = (85) 
The obvious solution of the last equation satisfying the boundary conditions (82) is: 

Xi,o{x) = Aq sm{^2Eix) 
v/2^a = m7r, m = l,2,... (86) 
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d=y 



L, 



b-d 



s=x 
B, 



Figure 8: The two bouncing mode bundles Bi and -B3 of the regular vertical skeleton in the 
rectangular billiards 



Coming back to the second of the equations (15) we can conclude that Xi, 1(2/1 ^) again is 
independent of y. 

Passing next to the third of the equations (15) and repeating arguments similar to those 
which led us to (84) we get the following equation for Xi,i{x)- 

^-^0^ + 2E,xiAx) + 2E2Xifl{x) = (87) 

with the solution: 

Xi,i(x) = Ai sin(v/2^x) + Bi cos(v/2^x) + ^2^^ cos(v/2^x) (88) 

The boundary conditions Xi,i(0) = Xi,i(fl) = enforce however Bi = E2 = 0. 
Using again (15) and the inductive arguments we come to the conclusion that xiiViX^X) 
is y-independent and the coefficients of its semiclassical series have the form: 



Xi,fe(x) =^fcsin(v/2^x), A; = 0,1,... (89) 
so is the form of Xi{x, A) itself, i.e. 



Xi(x, A) = ^(A) sin(V2^x) = ^(A) sin 



X 

mvr— 
a 



A{\) = J2 AkX'' (90) 

A:>0 
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Clearly, similar conclusion can be obtained for xsiu^ ^) which by (79) and for the mn-th 
energy level is equal to: 



X3,mn(y, X, A) = -(-l)"xi,m(x, A) = -(-l)M(A) sin (mTT^^ (91) 
Therefore coming back to (78) we get: 

*^^„(x, y, A) = 2iA(A) sin (^nvrQ sin (jmr^^ (92) 



which again is the result got in the previous section. 

The energy E is given however by the finite semiclassical series: 



E = y + ^ = ^ ^] I m,n = l,2,... (93) 



1 o El If /ra7r\2 /m7r\2 
2^ ^^ = 2[[jb) +iA^j ^ 

Let us stress the following main differences between the previous non-periodic case and 
the bouncing mode one despite the fact that in both the cases the results obtained are the 
same. 

1. The bouncing ball skeleton is represented in the phase space by a cylinder rather than 
by a closed torus; 

2. contrary to the non-periodic cases only one skeleton is sufficient in the bouncing mode 
case to get the whole spectrum of the energy; 

3. all terms of the semiclassical series expansion of GSWF exist (do not vanish) in the 
bouncing mode case, while only the zeroth one in the non-periodic one; 

4. the semiclassical series for the energy contains two first terms in the bouncing mode 
case and only zeroth non-vanishing term in the non-periodic one; and 

5. the JWKB approximation of the energy does not coincide with its global value. 

However similarly to the non-periodic case the bouncing ball solution is also regular and 
exact. 



5.3.2 Periodic skeletons different than the bouncing ball ones - singular SWF's 

Consider now skeletons which periodic rays do not bounce between the sides of the 
rectangle. The simplest such a case the fifth one in Fig. 7 is shown in Fig. 9. As in the non- 
periodic case there are again four bundles in the corresponding skeleton but contrary to the 
case mentioned only two rays (of four of them) belonging to two different bundles can meet 
at each point of the rectangle if this point docs not lie on SD's. 

We have to note also that the skeleton associated with this on Fig.9 differs from it by 
the opposite directions of rays, i.e. possible energy levels we get for these skeletons must be 
degenerated. 
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Figure 9: The singular skeleton defined by the corresponding pair Li and L2 of SD's 



The GSWF corresponding to the case looks as follows in different domains of the rectan- 



gles: 



oiApd3+iAps3 cos a 
^i\pdL+i\psL cos « 



X3(d3,53,A) + e^^K+^Ap4-n«^2(4,4,A) (x,y)Gl)3 (94) 
X3(d'3, ^3, A) + e^^P'^*+^^^^*--"x4(d4, ^4, A) (x, y) G 1)4 



, tana 



where the variables s^, are measured from the left ends of the corresponding sides A^, k = 
The Dirichlet boundary conditions on the respective sides of the rectangle are therefore: 



1,...,4. 



X4 



COS a 



,b — s tan a, A 







JXpa cosa^. 



COS a 



, s, A + X2(0, (a — s) tan a, A) = 



X3(0,s,A) + e^^f^^^""X2 (a -s) tan a, a) =0 

V cos a / 

f a — s 



V cos a ' 



, s, A + ^4(0, h — s tan a, A) = 



< s < a 



One can easily find from (95) that: 



Xi(0,s,A) = e2*^P(*'^'°"+'''=°'^")xi (^^,s,\\ 

\ cos Q / 



or 



g2i Ap(fe sin Q+a cos q) 



(95) 

(96) 
(97) 
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because xi{d,s,X) is periodic with the period 
the length of the rectangle diagonal. 



2(6 sin a + a cos a) = 2D where D is 



Therefore we get the following quantization condition for the zeroth energy term Eq = ^p^ 

XpD = mr, ra = l,2,... (98) 



Now we have to note that none of the bundle considered has a piece of its boundary 
common with any other one inside the rectangle. Seemingly as such could be considered 
the rectangle diagonals if the rays in the respective bundles were not run in the opposite 
directions. Therefore GSWF's defined in these bundles have to vanish on the respective 
diagonals of the rectangle so that we have to have: 



1,...,4 



Xk{d,0,X)=0, k 

0<d<D (99) 

But then from (95) we get also: 

XI (0, a, A) = X2(0, b, A) = X3(0, a, A) = X4(0, 6, A) = (100) 
Further using the propagation formula (15) for xi,i(2£),s) = xi,i(0)*) we get: 

X'iois) + 2Ei sin^ axi,o(s) = (101) 
which with the conditions (99)-(100) for xi('^, s, A) gives: 

Xi,o(s) = ^osin(v/2^sinas) (102) 
where Ei defines the second term of the semiclassical energy expansion with the condition: 

1 m^vr^ 



El 



m = 1,2, 3, ... 



(103) 



2 sin^ a ' 

Next repeating the procedure for the bouncing mode skeleton to the remaining terms 
Xi,ik(2D, s) = xi,fe(0, s), k = l,2, 3, we get for them: 

Xi,fc(s) = Ak sm{^/2E^ sin as) (104) 

so that 



Xi{d,s,X) = ^(A) sin(v2^ sin as) 



(105) 



ifc>0 

Therefore using (95) the final form of the SWF ^'"*(x, y) can be written as follows: 

' giAp„{xcosa+ysina) ('m7ri(x - y COt a) 



mn 



^iXp„(x cos a-y sin a) frmrKx + y cot a) 

^i\p„{x COS a+y sin a) g-^ (^^^^(a; - y COt a)) + 



^iXp„{-x COS a+ysina+2a COS a) (mTr^X + y COt a 
_giAp„(-a;cosa-2/sina) g^^^ {rmT^{x - ycoto)) + 



^i\pni-xcosa+ysina+2acosa) (m7r^{x + y COt a 
_^iXp„{-xcosa-ysina) (^m7ri(x - y COt « 
^i\pn(x cos a-y sin a) g^^ (mTT^ix + y COt 



{x,y) G 
{x,y) G D2 
{x,y) G 
{x, y) G £>4 



(106) 
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while the energy spectrum is: 



^ = 2^ + = 2)? + Wj ' - = 1' 2, 3, ... (107) 

By their construction the solutions (106) are all singular - their derivatives are discontin- 
uous on the rectangle diagonals. The corresponding energy spectrum also differs from the 
exact one as well as from the one of the regular SWF's obtained in the last two sections. 

It should be stressed however that the solutions (106) - (107) are allowed semiclassical 
solutions to the rectangle billiards eigenvalue problem, which corresponds physically to effects 
of the short wave limits. It is therefore of great importance whether one can detect in these 
limits resonant modes in the respective rectangular cavity corresponding to the SWF's (106) 
and the energy spectrum (107). Indeed such modes have been detected experimentally by 
Bogomolny et al [18] for the rectangle cavity with a barrier inside. This case of the billiards 
will be discussed in the next sections. 

If however such modes can be detected in the rectangle cavity then one can expect the 
corresponding GSWF's to have forms of standing waves rather than of the running ones as 
in (106). We can get such forms of GSWF's noticing that the spectrum (107) is obviously 
degenerate since in the case considered the associated skeleton is different from B. We 
can use therefore the corresponding running solutions for the skeleton B"^ to construct by 
superpositions the standing SWF's corresponding to the energy spectrum (107). The simplest 
two superpositions are: 

. / 1 . n\ I sm(Xpn(x COS a + y sina)) 

sm mTT^ia: — ucot a) X < , /, / ■ w — 

\ 0.^ " ^1 ± cos(Ap„(a: cos a + y sm ctj) 

. / 1/ n\ I sin(Ap„(xcosQ; — usina)) , , ^ ^ 

sm m7ri(a; + 2/cota) X <^ \/^ ) ^ . ii (x, w) e Di, L»4 

V "^1 cos(Ap„(xcosQ;-ysma)) ^ '^^ ' 

. / 1/ , n\ I sin(An„(x cos a + u sin a)) 

sm m7ri(x -j/cota) X <^ , /? / ■ \^ - 

V " 1 ^cos[Xpn(xcosa + ysmaj) 

. ( \, n\ I sin(Ap„(xcosa — -ysina — 2acosa)) / x -r^ 

sm m7ri(a; + ycota) X <^ / . ^^ (x,y) E D2, D3 

\ I — cos(Ap„(a;cos a — ysma — 2acosa)) ^ 

where the plus corresponds to the domains Di, D2 and the minus - to D3, D4. 

The above results on the singular SWF's can be generalized to arbitrary periodic skeletons 
in the rectangular billiards some of which are shown on Fig.6 and Fig.7. This can be done by 
noticing that an arbitrary periodic SD is defined by arbitrary two relatively prime numbers 
{p, q} so that a SD in a rectangle with the sides a and b shown in Fig. 9 starting from the 
vertex (0, 0) and being inclined by an angle a to the x-axis is defined by such two numbers 
as follows: 

tana=— (109) 
qa 

The above fact follows directly from the unfolded forms of periodic trajectories shown 
in Fig.6 and Fig.7 if one realizes that each of them has to finish on another vertex of the 



(108) 
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rectangle. It follows also that the set of all SD is countable but dense among all trajectories 
in the rectangle. 

Pairs {p,q} can appear in the following combinations {e, o}, {0,0} and {o,e} where e 
stands for "even" and o - for "odd". The respective SD's defined by these combinations 
finish their runs in the vertices (a,0), {a,b) and (0,6) correspondingly. 

A SD defined by a pair {p, q} bounces 1-times from each horizontal side of the rectangle 
and q — 1-times - from each of the vertical ones. If D denotes its global length measured from 
its starting vertex (0, 0) to one of its final ones just enumerated then D = pb sin a + qa cos a = 
ViqaJ^ + {pby^. 

If a SD is chosen, i.e. {p, q} are fixed, and it ends at one of the vertices just enumerated 
then the second SD which has to accompany the chosen one to built the skeleton starts and 
ends at the remaining two of these vertices. Note that a number of bundles in such a skeleton 
is then equal to 2p + 2q while their widths are equal to ^ sin ct = ^ cos a. 

The quantization formula (98) remains then valid for the case considered while (103) 
takes the form Ei = | ^ sirJa ' ~ 1)2,3,..., so that the formula (107) for the energy 
spectrum remains also unchanged. This formula corresponds to the spectrum of the exact 
standing wave functions in a rectangle with the sides D x ^sina. According to Bogomolny 
and Schmit [17] this rectangle can be considered as an unfolded skeleton so that each such 
a standing wave function defined on it should generate a corresponding semiclassical one 
defined on the skeleton by folding appropriately the rectangle mentioned to the real skeleton 
and interfering pieces of the standing wave function in crossed points of such a folding. 
However such a procedure to be correct still needs for the resulting SWF's to vanish on the 
rectangular billiards sides. It is seen that such a procedure though theoretically possible 
and correct is complicated enough to be replaced by the corresponding constructions of the 
SWF's in the real folded skeleton according to the rules formulated in sec. 3. 

Nevertheless the SWF's corresponding to the energy spectrum (107) are all singular having 
discontinuous derivatives on the lines separating two neighboring bundles. The lines are just 
the bundles boundaries on which the running SWF's defined in these bundles have to vanish. 

5.3.3 Regular SWF's in the periodic skeletons in rectangular billiards 

The singular SWF's found in the previous section provides us with the energy spectrum 
which is different from the regular one. However for a particular relations between the 
rectangle sides these solutions can become regular. This can happen if one assumes the 
following two additional conditions: 

\pb sin a = kir 
Xpa cos a = In 
A;,/ = 1,2,.... (110) 

where p is the global momentum of the billiards ball. 

While it is tedious to be checked the singular SWF's satisfying the conditions (110) become 
then regular obtaining the following form: 

.T.as I \ . (m-l \ . fm + k \ . . (m + l \ . (m-k \ /-.^-.^ 
= I — ^pvrx I sm ( — - — qiry j ± sm l——pTrx j sm ( — - — qiry j (111) 

valid in the whole rectangle area. 
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An immediate consequence of the conditions (110) is that they hmit the form of the 
rectangles for which the solution (111) can exist. Namely we have to have: 



tan2a=('^V = ^ k,l = 1,2,.... (112) 
\qaj I 

Let us make further an important note that despite its semiclassical origin the solution 
(111) is still exact being a linear combination of two solutions given by the formula (50). But 
since it is an eigenfunction of a given energy this combination means that both its terms have 
to be eigenfunctions of the same energy, i.e. we should have: 

(m - l)'^p^ (m + fc)2g2 _ (m + l)'^p^ [m - k^q^ 

2 ' U ~ l2 ' U (11^) 



a 



62 a? 62 



It is easy to check that this is the case if one takes into account the condition (^6)2/ = 
{qa)^k which follows from (112). 

Therefore if /cq, for which a satisfies (112) arc relatively prime then the remaining 
allowed pairs of k, I satisfying the condition (112) are of course of the form k = nk^, I = 
nlo, n = 1,2,..., and the energy spectrum formula (107) for such a rectangle takes the 
following final form: 

E = I^^ ^-^' + -'^^[f^^ + ^o)P\ m,n = 1,2,3,... (114) 
2A^ a'^ko 

Each energy level (114) is degenerate with the following base in their two dimensional 
degeneracy space: 



, , . {m-nlo \ . {m + nko 
^i,mn{x,y) = sm I pTTx I sm I qny 

^ , s . /m + nlo \ . [m- nko 

^2,mn{x, y) = sm ( piTX I Sm I 



m,n = 1,2,3,... (115) 

To conclude for the "periodic" skeleton considered in this section to get regular SWF's 
it is necessary for a rectangle to satisfy first the constrain (110). But if it happens then 
the corresponding SWF's coincide with the ones built on "generic" skeletons considered in 
sec. 4.1. An additional result of the considered periodic skeleton configurations is that the 
corresponding energy levels are degenerate. 

Let us stress however that if there are no such four natural numbers k, l,p, q by which for a 
given rectangle the condition (112) can be satisfied then none periodic orbit skeleton provides 
us with a possibility to construct on it regular SWF's except the bouncing ball skeletons. For 
these two cases the energy spectra obtained coincides with the one got from the "generic" 
skeleton calculations and their degeneracy disappears. 



5.4 Quantization of pseudointegrable systems - broken rectangles 

By a broken rectangle we mean the one which can be decomposed into a finite set of 
disjoint rectangles, see Fig. 10. If reintegrated it shows some number of rectangular bays and 
peninsulas. 
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Figure 10: An "arbitrary" broken rectangular billiards 



In fact the broken rectangles can serve as archetypes of pseudointegrable systems with 
an arbitrary genus. Since however we are interested in considering some special SWF's 
configurations related to classical periodic trajectories we shall limit ourselves to rather simple 
forms of the broken rectangles. The simplest one with a single peninsula (the L-shaped 
pseudointegrable billiards in terms of Kudrolli and Sridhar [2]) is shown in Fig. 11 and also 
in Fig. 12 and Fig. 13 with several skeleton configurations related to some SD. 

5.4.1 Regular SWF's in the broken rectangular billiards 

Skeletons in the broken rectangular billiards which are to provide us with the SWF's 
which would be the exact solutions to the corresponding eigenvalue problem have to be 
regular. Considering the billiards of Fig. 11 it is clearly seen that there are no such skeletons 
- each "generic" skeleton sketched on Fig. 11 has to have the bundles, denoted by 5. and 
6. in the figure, which can be composed with the bundles 2. and 4. respectively into two 
global but singular compound bundles. Therefore each "generic" skeleton has to be singular. 
Despite this the Lagrange surface form by these skeletons are closed and of genus 2. 

The first general conclusion which follows for the considered case of the broken rectangular 
billiards and the more so for the more complicated ones of Fig. 10 is that one cannot expect 
the obtained SWF's to be exact. 

To convince oneself of the correctness of the last conclusion we will consider the bouncing 
ball modes skeleton of Fig.l2C,D, instead of making a tedious calculations for the generic 
skeletons leading however to the same results. To this goal it is enough to match the corre- 
sponding GSWF's defined on the skeletons shown in Fig.l2C,D according to the conditions 
(32). Both the GSWF's have the form (92). Therefore the respective procedure leads us to 
the following quantization conditions for the energy Enm- 
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Figure 11: A single bay rectangular billiards with a "generic" global singular skeleton com- 
posed of six bundles. Its reduced form contains four global but singular compound bundles 
which form in the phase space the closed Lagrange surface of genus 2. 





Figure 12: A single bay rectangular billiards with singular skeletons built on periodic tra- 
jectories - the singular bouncing ball cases. Every skeleton forms the cylinder-like Lagrange 
surface in the phase space. 
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m,n = 1,2, . 



(116) 



where A^, Ay are the wave lengths of rays in the horizontal and vertical skeletons respectively 
shown in Fig.l2 and no, mo are the smallest integers satisfying Iga = uqc and kod = mob where 

^0,^0 Sire also integers. 

The respective SWF's are the following: 



where -D^r denotes the domain of the x, y-planc occupied by the broken rectangular of Fig. 11. 

The above GSWF is of course regular. Nevertheless due to the properties of the bouncing 
ball skeleton it is not exact despite the fact that it satisfies the Dirichlet boundary condi- 
tions as well as the SE. The reasons for its approximate character are the conditions (116) 
which cannot be satisfied if the corresponding length o, b, c, d of the broken rectangle are not 
commensurate by pairs. However if such a commensurability is satisfied by the sides of the 
broken rectangle then the solution (117) is exact. 

One can also easily realize that the last results can be generalized to any bouncing ball 
modes skeleton in the broken rectangular billiards of Fig. 10. A little bit surprising is that 
the semiclassical formulae (116) for the energy and (117) for the wave functions remain 
unchanged for any such a billiards while a number of conditions the wave lengths A^; and A^ 
have to satisfy filling the vertical and horizontal skeletons by integer numbers of their halves is 
increasing respectively to munbcrs of bays and peninsulas forming the sides od such billiards. 
In fact for the corresponding SWF's the half wave lengths ^A^- and ^Ay considered as the 
units of lengths on the respective horizontal and vertical sides of the broken rectangle have to 
measure these sides by integers. It means of course that these sides have to be commensurate 
so theoretically such a condition excludes SWF's for most the broken rectangular billiards. 
Practically however since incommensurability in fact does not exists in real measurements by 
experimental errors one can always tune the corresponding waves to the real dimensions of 
the broken rectangles. 

5.4.2 Singular SWF's in the broken rectangular billiards 

The SWF's (117) seem to be the unique regular ones which can be constructed in the 
broken rectangular billiards, i.e. any other SWF's should be singular. Examples of them 
corresponds to all the skeletons shown in Fig. 12 and Fig. 13. The skeletons of Fig. 12 define 
singular SWF's which are identical with the ones of the formulae (116) and (117) except that 
there are no relations between the rectangular sides a,b,c,d. These modes were observed 
experimentally by Kudrolli and Sridhar [2]. 



A sin ^ sin 




2^ = Asm{^7rx)sm{^y) {x,y) G D,r 



(117) 



{x,y) ^ Dbr 
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Figure 13: A single bay rectangular billiards with singular skeletons built on the shortest 
periodic trajectories different than the bouncing ball ones with the energy spectra of the 
respective GSWF's. Note that n, m = 1,2, for every spectrum. The red arrows show the 
outermost periodic orbits of the corresponding skeletons. The skeletons form the cylinder- like 
Lagrange surface in the phase space each. 
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2\-\a+4 h' {csina-k cos « Y I 




2A \a +4h (csin« — /)cos«) / 



d a 



2\^\{a — clf+b^ h^cos a 



Figure 14: A broken rectangular billiards with barriers and with possible superscar skeletons 
and with the energy spectra corresponding to the respective GSWF's. In all the above 
formulae n,m = 1,2, .... The corresponding Lagrange surfaces in the phase space are cylinder- 
like. 



Even more spectacular are skeletons built by periodic orbits different than the bouncing 
ball ones shown in Fig. 13. These are just the skeletons which provide us with the singu- 
lar SWF's with properties described by Bogomolny and Schmit [17] as superscars and was 
observed also experimentally by Kudrolli and Sridhar [2]. 

The singular SWF's corresponding to the broken rectangle billiards shown in Fig. 14 (upper 
figures) were observed by Bogomolny et al [18]. In fact the authors mentioned considered the 
limit of the billiards when d — c — >• (lower figures). They studied experimentally the high 
frequency modes in a microwave cavity [18] confirming the existence of the superscar modes 
predicted earlier by Bogomolny and Schmit [17]. 

For a completeness we shall give below the form of the singular SWF's for the skeleton A 
of Fig. 13 together with their degenerate energy spectrum. 



sin(^;rfit7^(^-ycot 

{x,y) GDi, Dq 
{x,y) eD2, Dr, 




sin(^(x cos 
lb cos(^(x cos 

/ n TT / 



^l,2;mn l-*^' IJ) — 

a + y sin < 



a)) 
a)) 



^ ^ ^ a + y sin 

sin( ^ (x cos a — y sin a) ) 
cos(^(xcosa — ysina)) 



sin( ^ (x cos a + y sin a) ) 
lb cos(^(x cos a + y sina)) 
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siii( c-g'cota (^-^cota)) X 

si"( c-7cota (^ + ^cota)) X 
[ {x,y) e D3, D4 



sin( ^ (x cos a + y sin a) ) 
lb cos (^(.T cos a + J/ sin a)) 

sin(^(xcosa — ysina — 2c cos a)) 
— cos(^(a;cosa — y sin a — 2c cos a)) 



tan a = - 
c 

\pD = nw 





nrn- 



1 



mn — nP 



El 7r2 



n 



+ 



m 



2(csinQ; — dcos a)"^ 
m,n = 1, 2, 3, ... 



2" ' A2 2A2 \^L>2 (csina-dcosa)2y' 

where D = + is the length of the diagonal shown in the Fig.l3A. 

5.5 The rational polygon billiards other than the rectangular ones 
triangle and the pentagon billiards 



(118) 



the 



In this and in the next section we will made a short review of the billiards systems which 
have been widely [8, 10, 18, 19, 20] considered both theoretically and experimentally having 
mainly in mind their skeleton description. 

5.5.1 The equilateral triangle billiards 

These billiards which dynamics is integrable have been considered in past very often [8, 23]. 
From the skeleton construction view point one can similarly to the rectangular case built the 
"generic" skeleton, see Fig. 15a, as well as skeleton generated by periodic orbits, Fig.l5b-d. 

In the "generic" case of Fig. 15a the corresponding skeletons consist of twenty four bundles. 
The reduced skeletons however have them already twelve and such a number counts a set 
of the BSWF's which as it is seen from Fig. 15a have to interfere to built the GSWF. The 
skeleton is regular and the result of such a superposition is well known [8, 23] so we do not 
perform a corresponding calculations which leads to the exact eigenfunctions and the energy 
levels as well of the corresponding eigenvalue problems. 

It has to be stressed however that since the associated skeletons B"^ coincide with B then 
the energy level degeneracy if happens has to be of different origins. 

Nevertheless we will report here the corresponding results for the periodic skeletons of 
Fig.l5b-c. The skeleton of Fig. 15b is defined by the SD composed of the three sides of 
the triangle so that the period of the orbit is equal to three. The skeleton contains only 
three bundles covering the whole triangle each so it is regular. To get GSWF it is enough to 
superpose in each point of the triangle only three BSWF defined on the bundles. It is obvious 
also that for the skeleton shown in Fig. 15b its associated one does not coincide with it having 
the same energy, i.e. the corresponding energy spectrum is degenerate. Taking therefore 
two independent superpositions of their GSWF we get the following regular semiclassical 
solutions for the case: 



^mni^^y) = sin sin (^-^mryj - sin ^rmr{x + VSy) 



sin 



mr X 
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Figure 15: Skeletons in the equilateral triangles: a) - a "generic" regular one forming a torus 
in the phase space, b) - defined by periodic orbits (thick arrows) forming a Mobius band in 
the phase space, c) - defined by periodic orbits with a cylinder-like Langrange surface, d) - 
with longer periodic orbits, e) - in the amputated corner triangle. The latter can be neither 
integrable nor pseudointegrable. 
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-m7r{x — VSy) 
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(119) 



and 

^'^^(x,?/) = cos I -m-KX 



niry + cos 



-rmr{x + VSy) 



sm 



rnr X 



V3 



cos 



-mvr X 



sm 



nir X + 



(120) 



with the energy spectrum: 



Emn = ^("^^ + 3"-^) 

m, n = 1, 2, ... 



(121) 



for both the solutions where m, n are even or odd simultaneously. 

Note that since both the above solutions are regular they are exact as well as their common 
energy spectra. 

In the case of the periodic skeleton of Fig. 15c defined by the periodic orbit of the length 
\/3 (the double height of the triangle) there are four bundles of the skeletons and the cor- 
responding four BSWF's which are singular since the skeleton is singular. Therefore the 
GSWF is defined locally in the domains D^, k = 1,2, 3, of Fig. 16c being however continuous 
in DiU D2U D3. It is the following: 
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sin \ 2rmT-^] sin(2n7ra;) (x,y) G Di 

sin ^2m7r-^ j sin(2n7ra;) — sin mir (^x + j sin mr (^x + \/3y) j {x, u) £ D2 

(x + V^y)] {x,y)eD3 



sm 



rmr 



with the energy spectrum: 



(122) 



nvr 



27r2 2 „ 2x 

m, n = 1, 2, ... 



(123) 



Of course ^'^^„(x, y) = if (a;, y) ^ L>i U L>2 U L>3. 

While the above ^^„(x, y) satisfies Schrodinger equation its derivatives are not continuous 
on the boundaries separating the domains D^, k = 1,2, 3. This is why it is only semiclassical 
approximation to the exact solutions given earlier. 

The spectrum (123) are naturally degenerate. However a source of this degeneracy is 
the symmetry of the equilateral triangle. Namely, there are solutions with the spectrum 
(123) which can be obtained by rewriting the solution (122) in the new coordinate systems 
obtained from the present one by moving it to the remaining two vertices and rotating it by 
ib^ respectively. 

These new solutions can of course interfere with (122) and with themselves so that trying 
to stimulate the corresponding state in a cavity certainly such superposed states will be 
generated rather than the "pure" state (122). 

To isolate however the state (122) it is enough to remove one of the three corners of the 
triangle as it is shown in Fig.lSe. 

The skeleton corresponding to the periodic orbit shown in Fig.l5d contains eighteen bun- 
dles none pair of which can be done a compound one. The length of the orbit is equal to 

D = 3\/7 while the wideness of each bundle is equal to w = Therefore the GSWF's 

for this case are singular and are linear combinations of the eighteen BSWF's. The energy 
spectrum is degenerate at least because the skeleton shown and its associated do not coincide. 
Its form can be obtained using the following general formula for the energy spectrum of the 
periodic skeletons: 
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1 2 El iWvr^ 
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where D is the length of the orbit period and w is the wideness of the skeleton stripe. 
Therefore for the case considered we get: 



En 



27r2 



[tx? + 2\t?\ 



63A2 

m, n = 1, 2, ... 



(125) 



with some possible relations between m and n which can follow from the corresponding 
possible relations between ^-factors of the BSWF's defined on the skeleton bundles and 
which can be established by detailed constructions of the corresponding GSWF's. 
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Figure 16: The pentagon billiards and its simplest whispering gallery superscar skeleton 
with five bundles. A typical periodic orbit is shown (the red straight line in the unfolded 
pentagon) the limiting form of which is the five-pointed star orbit (the SD green straight 
lines). It degenerates into the inscribed pentagon orbit (the blue straight line) with the half 
of the typical period equal to 10 cos ^ ~ 8.090. Gluing the end segments indicated by AB in 
the unfolded pentagon we get a Mobius band in the phase space. 



5.5.2 The pentagon billiards 

The pentagon form billiards were also the subject of intensive studies of both theoretical 
and experimental [19, 20]. In the latter case the corresponding pentagon cavities were made 
of some dielectric media. In the very high frequency region the corresponding electromagnetic 
waves form different modes among which the whispering gallery one of Fig. 16 was the most 
prominent. The other pentagon modes shown in the paper of Lebental et al [19] are more 
difficult for an identification in terms of the corresponding skeletons also because of different 
boundary conditions the authors wanted to consider. 

Nevertheless in our paper we would like to distinguish other pentagon modes of SWF's. 
We will not however consider in details a "generic" mode because of its complexity. Namely 
in such a case one can find simply by hand that there are twenty compound bundles com- 
posing the corresponding skeletons and consequently the same number of BSWF's which 
have to interfere to get the GSWF's satisfying the Dirichlet boundary conditions say. These 
compound bundles however are regular each so the corresponding skeletons are regular and 
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Figure 17: Another simple singular periodic skeleton in the pentagon billiards with the half 
of the period equal to 3 cos jq + 2 sin | ~ 4.029. Its Lagrange surface is cylinder-like. 



the GSWF's built on such skeletons will be regular and by that will be exact solutions to the 
corresponding SE with the exact energy spectrum despite the fact that these results can be 
obtained by the semiclassical approach. Therefore instead of the "generic" skeleton cases we 
consider the singular ones generated by periodic orbits. 

The simplest of such cases is the whispering gallery skeleton shown in Fig. 16. It is quite 
easy to write the corresponding singular GSWF. Nevertheless we limit ourselves to quote 
merely the corresponding result for the energy spectrum. Namely we have: 

2 2 \ 

m n \ 




+ 



25cos2| singly 

m,n = l,2, ... (126) 

where similarly to the triangle case the number m, n are simultaneously even or odd. 

The spectrum (126) is of course degenerate. Note also that in the pentagon (white) center 
the corresponding GSWF's vanish identically. 

Another case of the singular skeleton shown in Fig. 17 reminds the rectangular bouncing 
ball modes and the similar triangle modes of Fig. 15. The corresponding energy spectrum is: 



2A2 \^(3cosf^ + 2sinf)2 sin^ ^ j 

m,n = 1,2,... (127) 

This spectrum is of course degenerate - there are five independent solutions with the same 
spectrum. In the pentagon cavity all these solutions can be stimulated simultaneously. To 
isolate at least one of them it is enough to desymmetrize the pentagon into its forms shown 
for examples in Fig. 18. For the case a) of the figure the corresponding energy spectrum is: 

Emn = \ (^^^^ JL^r.r..JL\2 + ~2 



2A2 \ (2cotf^j + cosf^j) 



10 ' 10 ^ 10 



m,n = 1,2,... (128) 



while for the case b) it is: 



""""2X2 l^(2tanf^ + 3cos^)2 ^ sin^ 
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Figure 18: The deformed pentagon billiards with the singular periodic skeletons having the 
half periods equal to 2 cot jq + cos jq ^ 7.106 - the case a) and 2 tan + 3 cos jq ^ 3.503 - 
the case b) 



m,n = 1,2,... (129) 

6 Billiards with chaotic classical motions - superscars and pe- 
riodic orbits 

6.1 Singular SWF's in the chaotic polygon based billiards 

It follows from the previous section that the idea of the skeletons seems to be effective 
in solving some simple situations of quantum phenomena related semiclassically with bil- 
liards which shapes stimulate rather chaotic than regular (integrable or pseudointegrable) 
motions. An example of such cases is shown in Fig.lSe. Still more spectacular situations 
exist in billiards which can be obtained from the rectangular and the pentagonal ones by 
their deformations. Examples of such deformations and the superscar modes possible to be 
detected in such chaotic billiards are shown in Fig. 18 and Fig. 19. To describe analytically 
the superscar skeletons shown in these figures the methods of the previous sections can be 
applied directly. Note that the superscar mode corresponding to the Sinai billiards of Fig. 19 
was observed experimentally by Kudrolli and Sridhar [2] and by Sridhar and Heller [21] who 
studied the Sinai billiards also numerically (see also [5]). 

7 Summary and conclusions 

In this paper we have formulated very thoroughly the idea of skeletons in billiards and we 
have shown in great details how to construct on them semiclassical wave functions basing on 
the modified Maslov - Fedoriuk approach [14]. The modification mentioned utilizes rather 
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Figure 19: The deformed rectangular, the Sinai and the pentagon bilhards giving rise to 
chaotic motions with some possible superscar skeletons 



the corresponding configuration spaces than the phasespaces making use of the complex time 
to continue SWF's by caustics [22]. 

In general the skeletons in billiards play a role of the Arnold tori [15] in the integrable 
dynamical systems with smooth motions. 

In this paper we have considered mainly the rational polygon billiards beginning with the 
rectangular one to show that: 

• a huge variety of skeletons can be found in such billiards; 

• by the Keller - Rubinov [23] and the Kitchens - Berry [8] constructions global skeletons 
form in the phase space the closed Lagrangian surfaces of some genus (> 1); 

• the singular skeletons which are not global are represented in the phase space by open 
surfaces with boundaries. Most of the periodic skeletons belong to such a class and 
form in the phase space the cylinder- like or the Mob ius- like bands; 

• the SWF's which solve the energy eigenvalue problems in such billiards can be of two 
kinds - the regular and the singular ones - depending on the skeletons which they are 
constructed on; 

• the regular SWF's solve the energy eigenvalue problems exactly; 

• the singular SWF's solve the energy eigenvalue problems only approximately providing 
us with the superscar solutions of Bogomolny and Schmit [1, 2, 3, 17, 18, 19, 20]; 

• the superscar skeletons and SWF's can be found in many chaotic billiards which contain 
flat boundaries in favourable patterns [2, 3, 5, 21]. 

Some general conclusion which can be done by reassuming our results is that in addition 
to the saturation of the Gutzwiller formula [24] for the semiclassical Green functions the other 
role of the periodic orbits is their being of a local organizer of order in the integrable as well 
as in the pseudointegrable and the chaotic motions realized by the respective skeletons. 
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Appendix A 

In this appendix we are going to show, that the geometrical optics rule of reflections of 

rays off the billiards boundary is a consequence of demands of vanishing on the boundary of 
the linear combination (28) accompanied by the conditions (29) and (30). Namely consider 
the following superposition of SWF's: 

^l'{x,y;u,l;\) = ^l]-^{di, sv,u,l; X) + ^l^^{d2, S2;u,l; X) = 
J'l (di, si;u, /)e''^iH<ii+/;' cosa,,i(.';«,0<i.')4i^(di, si-u, I; A) + 

J^l id2, S2;u, l)e''''K'^'+fo' ^°^'^'=-2(^>'')'^^');^^22(d2, S2; u, I; X) (130) 

with 

r = [x,y] = rk,i{di,si;u,l) = ro(si) + di{si;u,l) = 

TCk,2{d2,S2-,U,l) = ro(.S2) + d2(.S2;n, /) 
rjk,i(ci, s; u, I) e Bk{u, I), Vk,2{d, s; u, I) e B'f.{u, I), Bk{u, I) / B',^{u, I) (131) 

i.e. the SWF's '^^\{di,si;u, I; X) and '^^\{d2, S2\ u, I; X) are defined respectively on the bun- 
dles Bk{u, I) and B'f,{u, I) with Dk{u, I) H D'f^{u, I) interfering in the crossing point [x, y] 
of two rays r^. i(di, si; n, /) and 2(t^2) ■52) belonging to the respective bundles. 
Therefore the condition for \l'^*(x,y, A) to vanish on Ai^{u,l) is: 

(0, S)e^'=^i /o cos akAs';u,l)ds'^a^^^^. ^ 
Jk,l (0' ^y''"'^ ^° cosak,2{s';u,l)ds' ^.2^^^. A) = 

ro{s) e Ak{u,l) (132) 
Because of the fc-dependence the last relation can be satisfied if and only if: 

rs rs 

a\ \ COS ak,i{s' ; u, l)ds' = a2 / cos ak,2{s' u, l)ds' , ro{s) E Ak{u,l) (133) 
Jo ' Jo ' 

It is easy to see however that there are only two solutions of the last condition: 

ak,i{s;u,l) = ak,2 {s;u,l) for ai = 02 
ak,i{s;u,l) = Tr - ak,2{s;u,l) for ai = -a2 

ro{s) e Ak{u,l) (134) 

The first solutions are however uninteresting identifying the bundles in a given segment 
and consequently leading to the solutions vanishing identically on Ak{u,l). 

Putting ak,i{s;u,l) = a{s;u,l) and ai = —a we get from the second solution and from 
(132): 

Xk^{s; u, I; A) = -Xk,i{s; u, I; A) = Xk{s; u, I; A) 

ro{s) e Ak{u,l) (135) 
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so that the combination (28) becomes: 

^l'{x, y; u, I; A) = si;u, I; A) - ^^^^5(^2, S2; u, l; A) = 

j;} (di, 5i; Z)e''^'=('^^+/o ' ^(di, ,1; n, Z; A) - 

ro(s)GAfe(n,0 (136) 

where Xkf^id, s; u, I; A), a = =b, are given by (27) with x^(0, s; u, I; A) = Xfc(*; '^^ h 

The last result shows that ^'^*(x, y; u, Z; A) vanishing on Ak{u,l) has to be represented 
semiclassically by a combination of at least two SWF's of opposite signatures and such that 
if ^((i, s; 7i, Z; A) is defined on the bundle Bk{u,l) then the second SWF "^'^2{d,s;u,l; X) 
has to be defined on the bundle B^{u, I). 

Appendix B 

It is shown in this appendix that S/^-^Uj, Ij) from the formula (38) is s-independent. To this 
end consider Fig.l on which a mapping hk{s;u,l) of the arc Ak{u,l) into an arc Af;i{u',l') is 
defined by the the bundle Bk{u, I). According to this mapping we have: 

ro{hk{s; u, I)) = ro(s) + D(s; u, I) 
ro(s) G Ak{u,l), ro{hkis;u,l)) G Ak'{u,l') (137) 

where D(s; u, I) is a vector linking the point ro(s) with the point ro(/ifc(s; u, I)) of the billiards 
boundary. 

Differentiating the equation (137) with respect to s we get: 

dhkis;u,l) djk{s;u,l) . 

cos P{hk{s;u, I)) = cos^(s) - D[s;u,l) sm'yk{s;u,l) + 

OS OS 

dD{s; u, I) 



cos jk{s;u,l)- 



ds 



dhkis;u,l) . • «/ \ , ,^^-fk{s■,u,l) 

smP[hk{s;u,l)) = sm^(s) + D[s;u,l) cos 7^(5; ii, /) + 

OS OS 

sin7.(.;«,/)^^^^ (138) 

where D(s; u, I) is the length of D(s; u, I). 

Making further the proper linear combinations of the last equations we have finally: 

dhk{s;u,l) , , dD{s;u,l) 
— cos ak'{hk{s;u, I );u ,l)= cosak{s;u,l) H — 

dhk{s;u,l) . n I ii\ ■ 1 n rx n^lk{s\u,l) 
sm. ak'{hk{s\u,ly,u ,l) = smak{s;u,l) - D{s;u,l) ^- (139) 

OS OS 

where we have taken into account the following relation between the angles involved: 

c^k'ihkis; u, l);u',l') + ak{s; u, I) = I3{hk{s; u, I)) - /3{s) + 2tt 

-fk{s;u,l) = /3{s) + ak{s;u,l) (140) 

which follows from Fig.l. 

The independence of s of 6k{uj,lj) follows now easily from the first of the relations (139). 
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